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ABSTRACT 


Performance  of  a vertical  array  for  estimating  multipath  signal 
parameters  and  for  detecting  multipath  Gaussian  signals  is  investigated 
in  this  report.  By  assuming  that  the  multipath  signal  directions  are 
known  and  the  number  of  array  elements  is  larger  than  the  number  of  paths, 
the  parameters  of  the  signal  from  each  individual  path  may  be  estimated 
so  that  the  multipath  signal  cancellation  problem  can  be  minimized.  The 
performance,  in  terms  of  signal-to-noise  ratio,  of  the  Maximum  Likelihood 
Estimate  (MLE)  for  the  signal  parameters  of  individual  paths  is  shown  to 
be  independent  of  the  relative  signal  phases  from  other  paths.  But,  the 
performance  does  depend  upon  the  angular  separation  of  multipath  signal 
directions  and  the  length  of  the  array.  In  the  case  where  some  knowledge 
about  the  signal  parameters  is  available,  the  Maximum  A Posteriori  Probability 
Estimate  (MAP)  is  also  considered.  The  linear  relationship  between  the  MLE 
and  MAP  estimates  is  derived.  It  is  interesting  to  point  out  that  the  MAP 
estimate  is  unique  even  if  the  number  of  array  elements  is  less  than  the 
number  of  signal  paths. 

The  optimum  array  processor,  for  detecting  multipath  Gaussian  signals 
with  known  arrival  angles,  has  a quadratic  structure.  In  a two-path  signal 
example,  the  optimal  array  has  two  beamformers  pointing  in  the  directions  of 
two  signal  paths  and  has  a term  which  cross-correlates  both  beamformers, 
even  when  the  signals  are  independent  from  path  to  path.  The  cross-correlated 
term  takes  care  of  the  leakage  of  signals  from  paths  other  than  the  beamforming 
direction  through  each  beamformer.  The  performance  evaluation  of  the  processor 
is  achieved  by  the  development  of  a general  method  called  "the  eigenvalue 
method".  This  method  may  be  used  to  find  the  analytical  expressions  for  the 
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probability  of  false  alarm  (Pp)  and  the  probability  of  detection  (P^O  for  array 
detectors  implementing  a quadratic  sufficient  statistic.  Several  examples  are 
presented  to  illustrate  the  use  of  the  eigenvalue  method. 

The  family  of  Receiver  Operating  Characteristic  (ROC)  curves  for  detecting 
multipath  Gaussian  signals  clearly  show  the  difference  between  the  performance 
for  detecting  multipath  Gaussian  signals  and  that  for  detecting  multipath 
known-waveform  signals.  The  ROC  curves  for  multipath  Gaussian  signals  cases 
are  not  only  a function  of  the  array  factor  and  the  multipath  structure,  but 
also  are  non-linear  on  normal-normal  probability  paper.  It  is  also  shown  that 
the  spatial  discriminating  capability  of  array  detectors  improves  performance 
at  high  signal -to-noise  ratios.  In  the  performance  comparison  between  optimal 
array  detectors  and  an  optimal  signal  array  element  detector,  we  found  the  array 
not  only  increases  the  output  signal -to-noise  ratio  by  k times,  where  k is  the 
number  of  array  elements,  but  also  improves  PQ  in  the  high  Pp  region. 

The  comparison  between  performance  of  an  optimal  array  and  a suboptimal 
array  is  also  carried  out  to  give  insight  regarding  the  trade-off  between 
processor  complexity  and  detection  performance.  In  the  two-path  signal  example, 
the  single  beamformer  array  which  points  a beam  in  one  of  the  two  signal 
directions  needs  10  elements  to  achieve  the  same  performance  as  an  optimal 
array  with  only  4 elements,  for  SNR=1  (signal-to-noise  ratio  per  element). 
Performance  of  a double  beamformer  array  which  has  two  beamformers  pointing 
in  the  directions  of  two  signal  paths  is  evaluated  for  comparison.  These 
results  show  that  the  performance  difference  between  the  double  beamformer 
and  optimal  arrays  is  noticeable  only  for  high  signal-to-noise  ratios  and  a 
small  array  aperture.  This  implies  that  the  cross  correlation  structure 
in  the  optimal  array  is  not  very  important  if  either  the  array  aperture  is 

(ii) 


large  enough  to  isolate  signal  paths  or  the  signal-to-noise  ratio  is  low. 

The  performance  comparisons  made  in  this  report  give  some  insight  regarding 
the  trade-off  between  processor  structure  complexity  and  detection  performance, 
which  are  major  factors  to  consider  in  implementing  array  detectors. 
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CHAPTER  I INTRODUCTION 

1.1  Underwater  Acoustic  Transmission 

The  acoustic  sound  propagation  in  an  ocean  channel  is  very  compli- 
cated. The  modes  of  sound  propagation  depend  on  the  sound  velocity 
profile  of  the  underwater  channel.  The  detailed  physics  about  the 
channel  characteristics  is  described  in  Tolstoy  and  Clay  [1]  and  Urick 
[2].  The  sound  velocity  in  the  ocean  depends  on  temperature,  pressure 
and  salinity.  Several  typical  velocity  profiles  are  shown  in  Fig.  1.1 
(from  Tolstoy  and  Clay  [1]).  In  general,  the  temperature  at  the  sea 
surface  is  higher  than  that  in  deep  water.  The  water  temperature  is  the 
factor  that  dominates  sound  speed  for  the  first  1,000  meters  in  depth. 
Beyond  that,  the  pressure  becomes  the  overriding  factor.  The  shape 
of  the  velocity  profile  determines  the  multipath  structure  of  the  sound 
channel.  Fig.  1.2  shows  the  sound  rays  associated  with  a particular 
velocity  profile.  The  sound  transducer  is  located  at  a depth  at  which 
the  velocity  is  minimum.  The  rays  are  calculated  using  a ray  tracing 
program.  By  this  method  we  may  predict  the  ray  incident  angle,  the 
transmission  loss  and  the  arrival  time.  Since  many  other  environmental 
factors  such  as  a random  surface,  a time-varying  sound  velocity  profile, 
and  a random  medium  are  not  considered  in  this  method,  the  results  ob- 
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tained  from  this  method  can  only  be  used  for  rough  estimates.  In  Fig.  1.2 
we  see  different  types  of  propagation.  Some  rays  which  do  not  touch  the 
surface  or  bottom  are  called  refracted- refracted  (RR)  rays.  Some  rays 


which  touch  the  surface  are  called  refracted-surface-reflected  (RSR) 
rays.  Some  rays  which  touch  both  the  surface  and  the  bottom  are  called 
surface-reflected-bottom-reflected  (SRBR)  rays.  The  multipath  signals 
received  at  the  receiver  depend  on  both  the  depth  and  the  range  of  the 
receiver  location.  This  result  has  been  presented  in  a paper  by 
Flanagan,  Weinberg  and  Clark  [3]. 

1.2  Review  of  Previous  Research 

Many  different  approaches  have  been  used  to  study  underwater  acoustic 
channel  characteristics  in  the  past  two  decades.  An  important  paper  by 
Steinberg  and  Birdsall  [4]  showed  that  the  phase  of  the  received  wave- 
form in  a fixed-system  in  the  straits  of  Florida  varied  less  than  100° 
during  intervals  of  1/2  - 1 hour  and  that  the  multipath  structure  of  the 
signal  had  stability  on  the  order  of  5 minutes.  Due  to  the  fact  that 
the  phase  is  much  more  stable  than  the  amplitude,  the  phase  estimates  are 
Fourier  transformed  to  find  the  relationship  between  the  internal  waves 
and  the  fluctuation  of  the  phase  spectrum  [5,6,7].  The  discontinuity  of 
the  phase  curve  which  is  associated  with  the  deep  amplitude  fading  is 
attributed  to  the  signal  cancellation  of  multipath  signals  [5].  The 
fluctuation  of  signal  energy  results  in  a nonstationary  confidence  inter- 
val for  the  phase  estimate.  Especially  during  a deep  fade,  the  phase 
information  is  lost.  A key  to  minimizing  this  signal  cancellation  is  to 
separate  the  multipaths  spatially.  This  will  be  discussed  in  later 
chapters. 

The  multipath  structure  of  the  underwater  acoustic  channel  was 
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received  signal  was  matched  to  the  sequence  to  find  the  multiple  time 
arrivals.  The  temporal  separation  of  different  arrivals  showed  the 
existence  of  multipaths  [4].  Dyer  [8]  examined  the  signal  fluctuations 
in  the  ocean  for  both  multipath  and  scattering  processes  and  showed  that 
the  multipath  propagation  dominates  the  fluctuations.  Adams  [9]  investi- 
gated the  fluctuations  of  the  transfer  function  of  a random  multipath 
channel.  Both  mean  and  variance  of  the  random  transfer  function  were 
calculated  for  his  multipath  model. 

In  order  to  obtain  spatial  information  about  the  channel,  linear 
arrays  have  been  used  to  process  underwater  acoustical  data  [11,  12,  13]. 

i 

Jobst  [10]  modelled  the  number  of  multipaths  as  a Poisson  distributed 
random  variable.  Using  this  model  the  spatial  coherence  and  temporal 
coherence  are  estimated  for  a horizontal  line  array.  It  is  shown  that 
spatial  coherence  decreases  with  increasing  angle  from  broad  side,  with 
increasing  frequency,  with  increasing  sound  speed,  and  with  increasing 
vertical  arrival  angle  at  the  receiver.  Williams  and  Battestin  [11]  used 
a vertical  beamformer  to  isolate  a single  RSR  (refracted-surface-reflected) 
path  and  to  investigate  the  phase  coherence  time  for  single  and  multipath 
signals.  Urick  [12]  showed  the  effect  of  multipaths  on  the  gain  of  an 

' 

additive  vertical  array.  The  phase  coherence,  in  terms  of  the  clipped 
correlation  coefficient,  falls  off  rapidly  with  horizontal  range  of  the 
source.  1‘owever,  no  amplitude  stability  of  a single  path  was  mentioned  in 
their  papers.  There  are  other  papers  that  investigate  theoretical  signal 
processing  techniques  for  arrays.  Hinich  [14]  used  match  filtering  concepts 
to  formulate  a maximum  likelihood  estimate  of  source  depth  using  a vertical 
array.  The  eigenfunctions  of  the  underwater  wave  guide  are  used  to  match 
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the  array  to  the  received  signal.  Bucker  [15]  used  a similar  method  to 
locate  sound  sources  in  shallow  water. 

1.3  Approach  of  This  Report 

The  relationship  between  the  error  in  the  phase  estimates  and  the 
signal-to-noise  ratio  is  presented  in  Chapter  II.  It  is  useful  to  have 
either  the  variance  of  the  estimate  or  a confidence  interval  along  with 
the  phase  estimate  so  that  we  may  judge  the  goodness  of  the  estimated 
phase.  We  will  show  that  multipath  signal  cancellation  can  contribute  to 
the  loss  of  signal-to-noise  ratio  during  a deep  fade  and  therefore  we  would 
like  to  investigate  the  signal  along  each  path.  To  distinguish  multipath 
signals  spatially,  a vertical  array  is  used  in  Chapters  III  and  IV.  The 
resolved  path  technique  is  introduced  in  Chapter  III  for  estimating  the 
signal  waveform  of  each  path.  Both  the  signal  amplitude  and  phase  are 
estimated  simultaneously.  The  signal  waveform  is  assumed  fixed  but  unknown. 
The  maximum  likelihood  estimate  of  the  signal  waveform  is  shown  to  be  an 
unbiased  linear  minimum  mean  square  error  estimate.  In  Chapter  IV,  each 
single-path  signal  is  assumed  Gaussian  with  known  variance.  Both  maximum 
a posteriori  probability  estimates  of  the  multipath  signals  and  the  likeli- 
hood ratio  test  are  formulated.  A beamforming  structure  in  the  optimal 
detector  will  be  pointed  out  in  Chapter  IV.  The  detection  performance  of 
optimal  array  detectors  is  evaluated  in  Chapter  V by  using  an  "eigenvalue 
method".  The  probability  density  function  of  the  sufficient  statistic  is 
shown  to  be  dependent  on  a set  of  eigenvalues.  The  general  form  of  the 
probability  density  function  is  derived  analytically.  The  probability  of 
detection  and  the  probability  of  false  alarm  are  also  obtained  in  closed 
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form.  The  eigenvalue  method  can  also  be  applied  to  evaluate  the  detection 
performance  of  suboptimal  detectors.  The  performances  of  suboptimal  pro- 
cessors, such  a single  beamformer  and  double  teamformers,  are  compared 
with  that  of  the  optimal  processor.  The  result  of  this  comparison  will 
clearly  show  the  trade-off  between  the  complexity  of  processor  structure 
and  the  detection  performance. 
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CHAPTER  II 


PHASE  ESTIMATION  PERFORMANCE  OF  A SINGLE  HYDROPHONE 

The  time  samples  of  signal  phase  estimates  have  been  used  for  investi 
gating  the  environmental  factors  in  the  underwater  acoustic  channel 
[4,5,6,7,11,12].  However,  the  performance  of  the  phase  estimate  usually 
is  not  presented  along  with  the  estimate.  Some  misleading  conclusions 
may  be  drawn  by  examining  only  the  phase  estimate,  especially  in  a deep 
fading  situation.  A typical  piece  of  underwater  data  is  presented  in 
Figure  2.1  (from  Steinberg  and  Birdsall  [4])  to  show  both  the  amplitude 
and  phase  of  the  received  waveform  simultaneously.  During  deep  fades, 
when  the  detector  output  is  less  than  -30dBvb,  the  phase  angles  have 
180°  or  90°  jumps.  The  deep  fades  may  be  attributed  to  multipath  signal 
cancellation.  Since  the  signal  in  that  experiment  is  a single  frequency 
sine  wave,  phasor  diagrams  are  shown  in  Figure  2.2.  A noise  free  signal 
phasor  diagram  shown  in  Figure  2.2(a)  demonstrates  signal  cancellation 
in  a four-path  channel.  If  signal  component  S^  rotates  150°,  then  the 
magnitude  of  the  resulting  signal  phase  changes  a lot.  The  resultant 
signal  phasor  in  Figure  2.2(b)  is  much  longer  than  that  in  Figure  2.2(a). 

The  fluctuation  of  signal  energy  does  affect  the  performance  of  the 
signal  phase  estimate.  Thomas  [17]  and  Cederquist  [18]  investigated  the 


error  probability  density  function  of  the  phase  estimate  as  a function 
of  signal-to-noise  ratio. 
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p(K»)  = exp 


cosij»  • exp 


erfc 


where 


erfc(X)  = 


VT  a 


COSi(/ 


and 


is  the  signal -to-noise  ratio. 


and  x is  the  error  of  the  phase  estimate. 

The  effect  of  signal-to-noise  ratio  on  the  phase  estimate  can  also  be 
demonstrated  graphically  in  Figure  2.3.  Suppose  there  are  two  signal 
phasors,  one  with  amplitude  a-j , the  other  with  amplitude  ag,  and  both 
are  corrupted  by  an  additive  noise  phasor  n^.  The  resulting  phase  estimate 
error  of  the  phasor  a-j  is  and  that  of  the  phasor  ag  is  In  Figure 

2.3,  we  can  see  that  ^ is  larger  than  As  shown  in  Figure  2.4,  the 
density  function  peaks  up  as  the  signal-to-noise  ratio  increases. 

Cederquist  [18]  also  calculated  the  limits  for  a 90%  confidence  level  [19]. 
The  curves  shown  in  Figure  2.5  are  from  Cederquist  [18].  We  can  see  that 
at  low  signal-to-noise  ratios  such  as  -lOdB,  the  90%  confidence  interval 
covers  ±160°.  This  implies  the  uncertainty  of  the  phase  estimate  is  very 
high  for  low  signal-to-noise  ratios. 

In  a single  hydrophone  case,  the  receiver  cannot  distinguish  signals 
with  different  incident  angles.  The  corresponding  signal-to-noise  ratio 
(SNR)  depends  on  the  relative  phases  at  the  hydrophone  of  the  different 
multipath  signals  as  shown  in  the  following  equation: 

Single  hydrophone: 

M 

Zv1*- 

m=l 


SNR  = 


where  am  is  the  amplitude  of  the  mth  path  signal,  and 
9m  is  the  phase  of  the  mth  path  signal. 

In  the  underwater  acoustic  channel,  the  variation  of  water  temperature 
due  to  tides  and  internal  waves  changes  the  acoustic  sound  transmission. 

The  phase  of  each  path  varies  widely  for  small  variations  in  the  sound 
velocity  profile  [1].  A detector  output  curve  shown  in  Figure  2.1  shows 
the  fluctuation  of  the  amplitude  of  signal  plus  noise.  If  we  assume  the 
output  signal-to-noise  ratio  fluctuates  as  much  as  the  detector  output 
shown  in  Figure  2.1,  then  it  is  meaningful  to  draw  confidence  intervals 
of  the  phase  estimate  along  with  the  phase  output  curve  shown  in  Figure 
2.1.  Clearly,  we  will  have  a nonstationary  confidence  interval.  It  is 
a reasonable  practice  to  consider  only  the  portion  of  the  phase  curve  in 
which  the  signal-to-noise  ratio  is  above  a certain  level.  Of  course, 
this  will  introduce  new  problems  in  dealing  with  data  discontinuities 
where  the  signal-to-noise  ratio  is  not  high  enough. 

A1  alternative  method  of  dealing  with  the  signal  cancellation  problem 
is  to  isolate  single  paths  by  using  a vertical  array.  In  Chapters  III  and 
IV,  methods  to  resolve  multipaths  spatially  are  investigated  with  the  hope 
that  the  resolved  paths  will  be  found  to  be  more  stable  and  that  other  infor- 
mation on  channel  characteristics  will  be  obtained. 


CHAPTER  III 

NON-RANDOM  PARAMETER  ESTIMATION  USING  A VERTICAL  ARRAY 


3.1  Problem  Geometry  and  Notations 

The  problem  geometry  of  a vertical  array  in  a multipath  channel 
is  shown  in  Figure  3.1.  The  time  delay  between  adjacent  elements 
associated  with  the  mth  path,  Tm>  depends  on  the  vertical  incident 
angle  $ . Since  the  incident  angles  of  multipath  signals  are  different, 
the  time  delays  of  multipath  signals  are  different  from  path  to  path. 

The  time  delay  associated  with  the  mth  path  can  be  expressed  as: 


D cos 


(3-1) 


where  D is  the  separation  of  adjacent  elements, 

C is  the  sound  velocity, 

a is  the  vertical  incident  angle  of  the  mth  ray. 


Throughout  this  report  the  incident  angles  of  the  multi  paths  are  assumed 

known.  In  a practical  situation,  these  incident  angles  may  be  obtained 
through  a ray  tracing  method. 

We  use  M to  denote  the  number  of  paths,  and  K to  denote  the  number 
of  array  elements.  We  consider  a single  frequency  sound  source  in  a 
multipath  channel.  The  received  data  is  Fourier  transformed  to  obtain  a 
phasor  of  a signal  frequency.  The  phasors  for  a time-delayed  waveform 

are  well-known  to  be  the  original  phasors  multiplied  by  a complex  phase 
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shift  e~JuT.  When  the  signal  is  present,  the  received  data  is  the  summa- 
tion of  multipath  signals  and  noise.  The  complex  vector  R represents  the 
phasors  received  by  array  elements  and  S represents  the  phasors  of  multi- 
path  signals. 

R = A S + n (3-2) 

where  II  = (Ri  ,R2,‘  • • ,R^) 

R^  : phasor  of  received  data  at  the  ith  element. 

sj  = (S^ .Sg.*** ,sM) 

Sm  : phasor  of  signal  from  the  mth  path. 
nT  = (n1  ,n2,*** ,nK) 

n^  : Gaussian  noise  received  at  the  ith  element. 

A = (Vq  »¥£»***  »Vf/|)  a K x M matrix 

\ - (1,  exp[ju>Tm],  exp[jwTm2],*** , exp[jVc(K-l )]) 
is  the  pointing  vector  of  the  mth  path. 

A four  path  noise-free  phasor  diagram  is  shown  in  Figure  3.2(a)  to 

illustrate  the  equation  above.  At  the  ith  element,  the  phasor  of  the  mth 

path,  S . is  equal  to  a eJ0m  where  am  is  the  amplitude  and  0 is  the 
r m ^ m m m 

phase  of  the  mth-path  signal.  The  signal  phasors  received  at  the  (i+l)th 
element  is  related  to  that  received  at  the  ith  element  by  the  following 
equation. 

S'  = S expCjoir] 
m m m 
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where  is  the  signal  phasor  from  the  mth  path  received  at  the  (i+l)th 
element  and  Sm  is  the  signal  phasor  from  the  mth  path  received  at  the 
ith  element.  The  time  delay  t contributes  the  phase  shift  <dt  . In 
Figure  3.2(b)  every  signal  phasor  in  Figure  3.2(a)  is  rotated  by  ajTm 
degrees.  The  resultant  noise-free  phasor  turns  out  to  be  larger 
than  R. . This  indicates  that  although  one  element  has  a potential  signal 
cancellation  problem,  other  elements  still  receive  some  signal  energy,  if 
the  vertical  angular  separation  between  paths  are  large  enough  or  the  array 
is  long  enough  to  resolve  paths. 


3.2  Maximum  Likelihood  Estimate  (MLE) 

When  the  signal  waveform  is  fixed  but  unknown,  the  maximum  likelihood 
estimate  can  be  used  to  estimate  the  received  multipath  signals  [201. 

Schweppe  [21]  developed  the  " decoup led-beam"  data  processor  for  estimating 
signal  waveforms  from  multiple-signal  sources.  The  maximum  likelihood 
estimate  (MLE)  used  in  his  paper  is  equivalent  to  the  Linear  Minimum  Mean 
Square  Error  Estimate  (LMSE).  The  technique  can  be  used  for  estimating  the 
multipath  signals.  The  MLE  is  derived  in  the  following  equations  by  maximizing 
the  conditional  likelihood  ratio,  a(R|$_). 


where 


A(R|S)  =*  exp[-(R-A  S)*2._1(R-A  S)  + R*^'1?.] 
a = E[n  n*]  * an2I 

a-1 3 -71 

an 


(3-3) 

(3-4) 


(3-5) 


S*  - [S, , S, 


• SJ 


Sm  is  the  phasor  of  mth  path  signal. 
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i\(R|S_)  Is  also  called  the  likelihood  function.  The  likelihood  function  is 
maximized  by  minimizing  the  first  term  of  the  exponent  in  equation  (3-3). 

Hence,  the  maximum  likelihood  estimate  of  the  signal  can  be  found  from  the 
folJowing  equation: 

R-AS=0  (3-6) 

If  the  dimensionality  K of  vector  R is  larger  or  equal  to  the  dimension- 
ality M of  the  vector  S,  then  we  may  solve  for  the  maximum  likelihood  esti- 
mate of  the  received  multipath  signal.  In  other  words,  if  the  number  of 
array  elements  is  no  less  than  the  number  of  multipaths,  then  it  is  possible 
to  resolve  every  path  and  to  find  values  of  the  MLE  for  signals  arriving 
along  each  path.  Otherwise,  if  the  number  of  array  elements  is  less  than 
the  number  of  paths,  then  the  MLE  estimate  is  indefinite,  (i.e.,  we  may 
have  an  infinite  number  of  solutions.)  In  general,  the  Maximum  Likelihood 
Estimate  of  the  received  multipath  signals,  S^,  can  be  expressed  as: 

iwLE  3 (AW*.  for  KiM  (3~7) 

★ 

Since  A A is  a K by  K matrix  and  the  rank  of  matrix  A is  not  greater  than 

v * * 

min(K,M)  the  inverse  matrix  of  A A does  not  exist  if  the  rank  of  A A is 
less  than  K.  This  implies  that  M cannot  be  less  than  K.  If  A is  a square 
matrix,  then  the  maximum  likelihood  estimate  for  the  received  multipath 
signals  is  quite  simple. 

SMI(r=A‘1R  for  K=M  (3-8) 
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The  mean  and  the  covariance  matrix  of  can  be  derived  for  K>M. 


= E[(A*A)'Vr] 

(3-9) 

3 (A*A)_1A*A  s 

= s 

(3-10) 

(3-11) 

covCS^leJ 

= e^^e"-^^mle‘^  ^ 

3 (A*A_)-1  ECd_  D.  j 

(3-12) 

* -1  2 

3 (A  A) 

(3-13) 

A K»2,  M*2  example  is  demonstrated  to  show  the  operation  of  the 
estimator  and  its  performance.  The  matrices  A and  A"1  are  shown  in  the 
following  equations. 


(3-14) 


(3-15) 


Since  the  matrix  A in  this  example  is  a square  matrix,  the  MLE  can  be 
expressed  as: 


1 


^-r 
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^MLE 


C3-16) 


1 

eJUT2  - eJa)Tl 


eja>T2  R]  - R2 
-ej<ijTl  R1  + R2 


(3-17) 


In  order  to  see  the  physical  meaning  of  the  matrix  operation, 
consider  a noise-free  case.  When  rv«0,  the  signal  phasor  received  at 
the  first  element,  R-| , is  equal  to  the  sum  of  the  phasors,  S-|  and  Sg. 
Similarly,  R2  is  equal  to  the  sum  of  the  phasors  S-j  and  S2,  where 
S^  a S.j  eJWTl  and  S£  3 S2  6^2.  The  noise-free  phasor  diagrams  of  §1 

a 

and  S2  are  shown  in  Figure  3.3.  The  phasor  components  of  R^  and  R2  are 
also  shown  in  Figure  3.  We  notice  that  the  resulting  phasor  R^  may 
be  very  small.  In  Figure  3.3(c).  the  phasor  diagram  of  the  maximum  likeli- 

A 

hood  estimate  of  the  signal  from  the  first  path,  S-p  shows  the  complete 
cancellation  of  the  signals  from  the  second  path.  Similarly,  the  phasor 

A 

diagram  of  S2  shows  the  complete  cancellation  of  the  signals  from  the 
first  path.  This  operation  Is  called  "infinite  sidelobe  rejection"  by 
Schweppe  [21] . 

The  mean  and  the  variance  of  the  estimate  are  usually  used  for 
estimation  performance  measurements.  It  is  quite  easy  to  show  that  the 
maximum  likelihood  estimate  is  an  unbiased  estimate. 


EtiMLEJ  = SC^ArVRJ  (3-18) 

= Et  CA*A)-1A*A.  S_  + CA*A)_1A*nJ  (3-19) 

3 & (3-20) 

A A 

The  variances  of  S-j  and  S2  are  also  derived  in  the  following  equations. 
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Varf  S-j  ] - E[(S1  - S ,)*]  (3-21) 

- 1 io>T  1 W 1^  [Var(R,)  Var(R2)1 

eJanr2  - eJ(JTl  2 1 2 (3-22) 


2<t 


| eJU,T2  - eJtlJTl  | 2 

Similarly,  we  may  show  that: 

Var[S2]  = ET(S2  - S2)2] 
s Varf  §•,] 


(3-23) 


(3-24) 

(3-25) 


From  Chapter  II,  we  know  that  the  phase  angle  estimate  depends  on  the 
signal-to-noise  ratio.  It  is  interesting  to  compare  the  performance  of 
the  signal  phase  estimate  of  the  resolved  path  method  with  that  of  the 
single  hydrophone  method.  The  expressions  for  the  signal-to-noise  ratios 
of  these  two  methods  are  shown  below. 


Resolved  path  (MLE)  : SNR 


ejwT2  . ejaJT1 


2cj 


for  i=l ,2 


(3-26) 


Single  hydrophone:  SNR  = 


2 


(3-27) 


a.,eJ0l  + a2eJ02 

where  a^  is  the  signal  amplitude  of  the  ith  signal  path, 
and  9.j  is  the  signal  phase  of  the  ith  signal  path. 

In  the  single  hydrophone  case,  the  signal-to-noise  ratio  goes  to  zero 
even  for  multipath  signals  with  large  amplitudes  if  the  signal  phases,  e1 
and  $2»  are  180°  out  of  phase  and  the  signal  amplitudes  are  identical. 

This  signal-to-noise  degradation  may  be  minimized  by  employing  the  spatial 
discriminating  capability  of  a vertical  array.  As  shown  in  equation  (3-26) 
the  signal-to-noise  ratio  of  the  resolved  path  maximum  likelihood  estimate 
is  a function  of  time  delays  t-|  and  t2>  where  the  time  delays  are  proportional 
to  the  array  element  spacing.  Since  both  and  t2  depend  on  the  separation 

of  array  elements,  proper  control  of  the  separation  of  the  array  elements  will 

2 2 

increase  the  signal-to-noise  ratio  to  its  maximum  at  2 aT/c^  . For  a 

2 2 

K-element  array  the  maximum  signal-to-noise  ratio  may  reach  K a./an  . From 
the  discussion  above  we  observe  that  the  signal-to-noise  ratio  of  the  resolved 
path  estimate  depends  on  the  spatial  resolving  power  of  the  array  and  the 
signal-to-noise  ratio  can  be  Increased  by  increasing  the  array  aperture.  If 
the  "noise  like"  fluctuation  of  the  signal  amplitude,  as  shown  in  the  single 
hydrophone  case  [4],  is  truly  due  to  multipath  signal  cancellation,  then 
the  resolved  path  signal  estimate  should  be  stable  enough  to  provide  additional 
signal  amplitude  information  for  studying  the  channel  characteristics. 

On  the  other  hand,  equation  (3-27)  shows  that  the  signal-to-noise  ratio 
of  the  single  hydrophone  estimate  depends  highly  on  the  relative  phases  of 
signals  from  different  paths.  As  shown  in  Tolstoy  and  Clay  [1],  the  variation 
of  signal  phase  due  to  the  variation  of  the  channel  temperature  profile  is 


different  from  path  to  path.  This  implies  that  the  summed  signal  phase 
may  fluctuate  significantly.  Due  to  the  potential  signal  cancellation, 
the  signal-to-noise  ratio  and  the  confidence  interval  of  the  phase  estimate 
may  vary  dramatically  in  the  single  hydrophone  case.  Since  the  fluctuation 
of  x ^ and  tg,  which  depends  on  the  incident  angles  of  planewave  arrivals 
are  relative  stable  compared  to  that  of  summed  signal  phases,  we  may  con- 
clude that  the  resolved  path  method  should  provide  additional  information 
for  studying  both  the  signal  amplitude  and  phase. 


CHAPTER  IV 


ESTIMATION  AND  DETECTION  OF  RANDOM  MULTIPATH  SIGNALS  USING  A VERTICAL  ARRAY 


The  maximum  a posteriori  probability  estimate  (MAP)  for  each  indivi- 
dual signal  path  is  presented  in  this  chapter.  The  relationship  between 
the  maximum  likelihood  estimate  and  the  maximum  posteriori  probability 
estimate  is  also  investigated.  The  optimal  array  signal  detector  which 
uses  the  likelihood  ratio  test  will  be  shown  to  have  a beamformer  structure 
and  cross-correlator  terms.  Only  the  structure  of  the  detector  will  be 
discussed  in  this  chapter.  The  performance  of  the  detector  will  be  presented 
in  Chapter  V. 


4.1  Maximum  A Posteriori  Probability  Estimate  (MAP) 

Suppose  the  multipath  signals  are  uncorrelated  zero-mean  Gaussian 
random  processes.  Then  we  may  find  the  maximum  a posteriori  probability 
estimate  (MAP)  for  the  random  multipath  signals.  Assuming  we  know  the 
covariance  matrix  of  the  received  multipath  signal  vector  S_,  the  MAP  estimate 
can  be  obtained  by  maximizing  the  probability  of  given  R.  The  expression 
for  p(S|R)  is  shown  below. 

ptsjfi)  ■ ttrj1-  <4-') 

■ <4-2) 


where 


p(R|S)p(S)  = C]  exp[-(R-A  S)\](R-A  i)  - iO^S] 

and  Qs»  ^ are  the  covariance  matrices  of  received 

multipath  signals  and  received  noises  respectively. 
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Schweppe  [21]  has  derived  the  maximum  a posteriori  probability  estimate 
for  Gaussian  signals  in  uncorrelated  Gaussian  noise.  The  maximum  a_ 
posteriori  probability  estimate  of  the  signal,  S^p.  is  expressed  in 
equation  (4-4).  Unlike  the  maximum  likelihood  estimate,  S^p  does  not 
have  a K > M restriction. 

= [A.  A + ] A B.  K = M (4-4) 

Although  the  number  of  array  elements,  K,  is  assumed  no  less  than  the 

number  of  paths,  M,  in  this  chapter,  the  expression  above  is  valid 

even  for  K < M.  We  recall  that  in  the  expression  for  the  maximum  likeli- 

★ 

hood  estimate,  the  matrix  A A cannot  be  inverted  if  K < M.  But,  in  the 

Sr  O _ 1 

expression  for  the  MAP  estimate  the  inverse  of  the  matrix  [A  A + 2^  ] 

do*»c  exist  even  for  K < M,  because  the  rank  of  the  matrix  is  equal  to 
max(K,M).  This  is  due  to  the  introduction  of  a_  priori  knowledge, 
into  the  processor.  So,  even  when  we  do  not  have  enough  array  elements  to 
resolve  every  path,  we  can  combine  the  a priori  knowledge  with  the  observed 
data  to  do  a MAP  estimate  for  the  Gaussian  signal  for  each  path. 

1-v.r  the  cases  where  K >_  M,  there  is  a relationship  between  the  MLE  and  the 
MAP  estimator.  The  expression  for  S^Ap  can  be  rewritten  as: 


S^p  » Cl  + an2(A*A)‘1gsr1  (AW!!  for  K > M (4-5) 
= CL+  °n2(A*A  )'12s]'%LE  (4-6) 


mi  ph^i  i ■ 1 1 fi warn 


T 


— 
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This  expression  shows  clearly  how  to  bring  in  the  a priori  knowledge  into 
the  MAP  estimator.  Since  all  matrix  operations  are  linear,  we  may  derive 

A ^ 

^MAP  from  ^MLE  by  a linear  transformation  for  K>M. 


4.2  Structure  of  the  Optimal  Array  Processor  for  Detecting  Multipath 
Gaussian  Signals  in  Gaussian  Noises 

The  optimal  array  detector  can  be  formulated  by  implementing  the 

likelihood  ratio  A(R)  which  is  defined  as  A(R)  = r 201 . The 

P(R|HQ) 

probability  of  R,  which  is  a complex  random  variable  given  the  signal  is 
present,  can  be  expressed  as  [22]: 


p(R|H1)  = Tr"K{det  expt-jfg^RJ  (4-7) 

where  * E[ R * A EfS  SJA*  + E[n  nj  (4-8) 

Similarly,  the  probability  of  R given  the  signal  is  absent  is: 

p(R|HQ)  = ir'K{det  Qy}'1  expC-jM^RJ  (4-9) 

where  Qg  « E[R  R* | HQ]  » E[n  nj  (4-10) 


The  likelihood  ratio  of  jl  is  the  ratio  between  equations  (4-7)  and  (4-9). 
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* -1  * _i 

n 1 n i n n » » 


W 3 dir^"  exrt  -a  0.!  R + R Qo  RJ 


(4-11) 


A sufficient  statistic,  y,  which  is  a monotonic  function  of  A(R),  is 
shown  in  the  next  equation. 


y 3 MQo'1  • v1]  i 


(4-12) 


To  simplify  the  detector  structure,  we  assume  that  the  multipath 
random  signals  are  independent  from  path  to  path. 


Then, 


2,  ■ Qq  + E[  A S S A ] 

* Qo  + - 


(4-13) 

(4-14) 


*%  + — cr 


O ■ 


(4-15) 


5o  + S gL  \ 


“i  s_  — m — m 
m*l  m 


(4-16) 
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where  A 3 [ Vf  .Vg,***,^  ] 

and  <j^  is  the  variance  of  signal  from  the  mth  path. 


(4-17) 


From  the  structure  of  as  shown  in  equation  (4-16),  we  may  apply  the 
well-known  matrix  inversion  lermia  M times  to  find  \ The  lenina  [23] 
is  stated  in  the  following  equation. 


11  1 -1  * -1 

[B  + u u*l  3 B S— T — B u u B 

“ (1  + u B’*u)  ~ " 


(4-18) 


where  u is  a vector  and  £ is  a square  matrix. 

A M*2  example  is  formulated  to  give  more  insight  into  the  detector 
structure.  In  t*e  example  we  assume  that  the  noises  are  independent  from 
element  to  element,  so 


a,  - 4 i * ^ \ W 


(4-19) 


We  let  the  first  two  terms  in  the  equation  shown  above  be  B 1 , such  that 
2 


B.  3 o?  I + a\  v,v,  . Then, 
— 1 n — s-j  —i— i 


2/4 

-12  2-.11  ,V» 

8.1  " i + \ 1 ' J -~J~3 

1 n s 


*1*1 


(4-20) 


1 + K 


1 

3T 


L 


After  some  manipulations,  we  find: 
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★ 

We  may  call  v,  R a beamformer  pointing  in  the  direction  of  the  first  path. 

★ 

Similarly,  Vg  R.  1S  a beamformer  pointing  in  the  direction  of  the  second 

"tt  "if  'tr  O 

path.  The  term  (R  v^)(v^  R)  = |v^  R|  can  be  physically  realized  as  a 
beamformer  pointing  in  the  direction  of  v^  followed  by  an  envelope 
detector.  The  structure  of  the  optimal  array  detector  is  shown  in 
Figure  4.1  as  a combination  of  beamformers,  envelope  detectors  and  cross 
correlators.  The  first  and  second  terms  in  equation  (4-30)  measure 
the  output  power  of  beamformers.  The  third  term  is  a correction  term 
which  correlates  the  output  of  two  beamformers.  The  effect  of  the  third 
term  on  the  detection  performance  will  be  discussed  in  Chapter  V. 

Although  we  only  dealt  with  a single  frequency  case  in  this  part,  it 
is  not  difficult  to  extend  this  result  to  broadband  signals.  The  difference 
between  the  single  frequency  and  broadband  signals  cases  is  that  one  uses 
time  delays  instead  of  phase  shifts  for  the  broadband  signals. 
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CHAPTER  V 

PERFORMANCE  OF  DETECTING  MULTIPATH  GAUSSIAN  SIGNALS  IN  GAUSSIAN 
NOISE  WITH  A LINEAR  ARRAY 

A general  method  called  the  "Eigenvalue  Method"  for  evaluating  per- 
formance of  array  detectors  will  be  introduced  in  this  chapter.  The 
method  is  general  enough  to  cover  all  detectors  which  employ  the  quadratic 
sufficient  statistic  (i.e.  y=R  Q_  R , where  ^ is  a hermitian  matrix.)  The 
performance  of  the  optimal  detector  will  be  evaluated  and  compared  with 
that  of  conventional  suboptimal  detectors  such  as  the  single  beamformer 
and  multi beamformer. 

5.1  The  "Eigenvalue  Method" 

Suppose  the  array  signal  detector  has  a quadratic  sufficient  statis- 
* , * 

tic  y = R £ R where  £ is  hermitian  (i.e.  £ = (£).  Then,  we  may  find  the 
characteristic  function  of  y using  the  following  equations. 

♦/1«>  = Ey[exp(i$y)]  (5-1) 

= ER(exp[  i<tR*Q_  RJ } (5-2) 

exp(i$R*(i  R)]p  (R)d  R (5-4) 

where  p(R)  = ir“k{det  V}"^exp[ -R  V-^  RJ  , (5-4) 

V = cov[R]  , 


and 


<t»y(i$)  is  the  characteristic  function  of  y. 
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Then,  #y(U)  n"k{det  V)‘1exp[-JlV,R  + i*R*£  R]d  R (5-5) 


, f i r'K(det  VJ^expf-RV1^  - i*V  gjRJdR  (5-6) 

» n 


{detfl  - i*  V 0]  >“ 


(5-7) 


To  simplify  the  notation,  we  let  G 3 V £,  where  G is  a square  matrix  with 
dimension  K.  Any  square  matrix  G is  similar  to  either  a diagonal  matrix 
or  a triangular  matrix.  We  may  find  matrices  £ ^ and  £ such  that 
£“^G  £ 3 £ where  £ is  either  a diagonal  matrix  or  triangular  matrix  with 
all  the  eigenvalues  of  G on  the  main  diagonal.  Then,  det[£  - i$G]  can 
be  expressed  as: 


detfl  - i$  £]  3 det{£"  }det[£  - i<t  G]det{£} 
3 detf  I - i<5  £_1£  £] 

3 detf  I_  - i$  £) 


(5-8) 

(5-9) 


(5-10) 


0 - i*Xa) 


(5-11) 


where  X^  is  the  eigenvalue  of  matrix  G, 
Hence,  the  characteristic  function  of  y becomes 


V1*)  s TT  (1  ‘ i*x*)~1 

l- 1 


(5-12) 


1 
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The  density  function  of  y can  be  obtained  directly  by  taking  the  inverse 
transform  of  the  characteristic  function, $y(i $)  of  y.  Dyer  [26]  used 
this  technique  to  derive  the  expressions  for  the  density  functions  of  y as 


(1)  All  X^s  are  equal  to  X. 


p(y) 


yK_1exp(-y/X) 


for  y>0 
for  y<0 


(2)  Each  x has  a distinct  value 

)C 


p(y) 


2 X."1  h "1  exp(-y/X.) 

l»l  1 

0 

0 


for  y>0 
for  y<0 


(5-13) 


(5-14) 


where  h^  = 


(5-15) 


(3)  m groups  of  eigenvalues,  each  group  has  identical  eigenvalues. 


exp(-i  k 


n<kt'V  1 


k 3 -iX_ 


(5-16) 


The  above  equations  are  derived  using  calculus  of  residues.  The  first 
two  equations  are  very  useful,  but  the  last  equation  does  not  give  us 
any  idea  of  the  functional  form  of  p(y).  Hence,  another  approach  is 
taken  in  this  report  as  will  now  be  shown. 

We  recognize  that  the  probability  density  function  corresponding  to 
the  characteristic  function  (1  - itX^)"^  is  exponentially  distributed. 

In  other  words,  if  4>x ( 1 <&)  * (1  - itX^)-1  then. 


p(x)  A r-  exp(-  x/X.) 


for  x>0 


for  x<0 


(5-17) 


We  can  view  the  random  variable  y as  a sum  of  independent  random  variables, 

x.  , because  the  characteristic  function  of  y can  be  factored  as 

4>  (i$)  = (i<&)*<f>„  (i<t) •••$„  (i4>).  Since  the  probability  of  x_.  is  expo- 

y x1  x2  xk  1 

nentially  distributed  as  described  in  equation  (5-17),  the  probability 


density  function  of  y is  the  convolution  of  exponentially  distributed  proba- 
bility densities.  The  resultant  function  p(y) , obtained  from  the  convolution,  will 
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take  one  of  the  following  three  functional  forms. 


(1)  All  X.  are  equal  to  X. 


1 K-l 


1W 


exp(-y/X)  for  y>0 


for  y<0 


(2)  Each  X^  has  a distinct  value. 


(5-18) 


P(y)  4 2Z  exp(-y/X^) 


for  y>0 


for  y<0 


(5-19) 


where  3 — 


- Xj) 


(3)  m groups  of  eigenvalues,  each  group  has  M identical  eigenvalues. 


p(y).|t  z\r/(M^'r  Wv]  f0r  y-  («) 

) t»l  L r*0  ri 


for  y<0 


where  is  the  number  of  elements  in  the  4th  group. 
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m 

E M*  a K> 

»-i 

and  d depends  on  the  eigenvalues,  but  is  not  a function  of  y.  From 
rt 

equations  (5-18),  (5-19),  and  (5-20),  we  can  see  that  p(y)  is  a finite 
sum  of  functions  with  form  ymea^.  This  functional  form  will  give  us 
closed  form  expressions  for  the  probability  of  detection  (Pp)  and  the 
probability  of  false  alarm  (Pp).  By  definition,  PQ  and  Pp  can  be 
written  as: 


P 


D 


(5-21) 


Pp  - f P Cy I Hq) dy  (5-22) 

n 


If  we  perform  these  integrals,  we  will  have  expressions 
which  can  be  found  in  an  integral  table  as: 


vm  ay 

J e 


dy. 


n 


m-r 


m _m-r 

ean  E C-D^1  1111 

r*0  (m-r)l  ar 


(5-23) 
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Eventually,  we  will  not  have  any  integral  left  in  the  expressions  for  Pg 
and  Pp.  In  the  next  two  sections,  we  will  find  the  eigenvalues  of  G along 
with  the  expressions  for  Pg  and  Pp  analytically,  for  some  special  cases.  In 
general,  the  eigenvalues  of  the  G matrix  can  be  calculated  using  numerical 
methods.  There  are  computer  programs  developed  at  Arqonne  National 
Laboratory  and  currently  available  at  many  computing  centers  [24] 
for  solving  eigenvalues  of  any  complex  matrix.  After  solving  for  the  eigen- 
values of  matrices  Gg  and  G^ , equations  (5-18),  (5-19),  (5-20),  (5-21), 
and  (5-22)  can  be  used  to  calculate  Pg  and  Pp.  This  method  which  is  called 
the  "Eigenvalue  Method"  can  be  used  to  calculate  the  exact  values  of  Pg  and 
Pp  at  any  point  on  the  ROC  curve  with  the  aid  of  a computer. 


5.2  Optimal  Performance  Evaluation  for  Array  Detectors 

In  this  section,  we  will  investigate  some  special  properties  of  the 
eigenvalue  \ ^ in  the  optimal  detector  cases.  We  will  also  derive  eigen- 
values analytically  for  some  special  cases.  Performances  of  the  detectors 
are  compared  to  see  the  effect  of  adding  hydrophones  to  the  array  receivers. 

First,  we  would  like  to  derive  a general  formula  for  performance 
evaluation.  From  Chapter  IV,  we  saw  that  a sufficient  statistic  for 
likelihood  ratio,  y,  can  be  expressed  as: 

y a - if' ' )R  (5-24) 


Qg  iS  COV[  R | Hg] 

9,-j  is  covf  R|H.|] , 


and 
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= Ey[exp(i<&y)] 


= ER{exp[i<fR*(^1-£j‘1)Rj} 


= rr“K{det 

= / ir~K(det  V}-1exp[ -R*V-1n  - iW^1  - g^1 )}R] dR 


V}_1expf  -rV1!*  + UR^gjj1  - g.^)RjdR 


3 {det(I  - i*V(g^  - 5.J1)]}'1 


where  V * cov[RJ 


(5-25) 

(5-26) 

(5-27) 

(5-28) 

(5-29) 


We  let  G = V (Q^1  - g^)  to  simplify  the  notation.  The  matrix  G under 
both  HQ  and  can  be  further  expressed  as: 


Hq:  cov[R|H0]  = ^ 

(5-30) 

% = QoQo1  "2.1^ 

C5-31) 

(5-32) 

u 

I-1 

i 

J? 

1 

(5-33) 

h,:  covcaiH,]  ■ a. 

’(5-34) 

Sl  ■ Si  (si1  - si’) 

(5-35) 

- a,  aS1  - i 

(5-36) 

mmmm 
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From  equations  (5-33)  and  (5-36),  we  may  find  the  relationship  between 
the  eigenvalues  of  G and  that  of  as: 


x 


= 1 


1 

* -1 


V0 


-1 


- 1 


(5-37) 


(5-38) 


These  two  equations  show  that  we  need  only  to  solve  for  the  eigenvalues  of 

2q  Qo"1  and  then  use  equations  (5-37)  and  (5-38)  to  calculate  Xg  and 

X-  . Furthermore,  we  may  derive  the  relationship  between  \r  and  x„  . 

1 G0  G1 


(5-39) 


(5-40) 


G 


(5-41) 


(5-42) 


This  implies  that  once  we  know  Xr  we  can  find  \r  very  easily. 

b0  G1 


35 


There  are  special  cases  in  which  analytical  expressions  for  the 

eigenvalues  of  G can  be  obtained.  These  special  cases  will  be  shown  as 

2 

examples.  To  simplify  the  derivation,  we  further  assume  that  (Jq  3 an  I_, 
which  means  the  noises  at  the  receiving  array  elements  are  independent 
from  element  to  element  and  have  equal  power. 

Example  1 : Single  Path  Signal 

If  we  have  only  one  path,  then  the  covariance  matrix  ^ can  be  written  as: 

2j  3 °n  I + (5-43) 

The  matrices  Gg  and  G^  become: 

So  - 1 - 4 a,*'  <5-44> 

■ i -<£<<£  i ♦ »!,  ii  O'1  <5-45> 

a,  ■ 1 4 - 1 (5-«) 


■4  S,  - I (5-«) 

a 


’ 
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ka 


*1 


2 2 
Ka/  + o‘ 

S1  n 


(5-54) 


From  equation  (5-49),  we  observed  that  the  rank  of  is  one.  By  using 
equations  (5-40)  and  (5-42),  we  know  that  a zero  eigenvalue  of  either  G^ 
or  Oj  will  imply  a zero  eigenvalue  of  the  other.  So,  the  rank  of  Gq  is 
also  one.  Hence,  the  density  functions  of  y under  both  hypotheses  are 
exponentially  distributed  as: 


p(y|H0) 


p(y  I H-, ) 


~ exp(-y/XG  ) for  y>0  (5-55) 


for  y<0 


~ exp(-y/x  ) for  y>0 


for  y<0 


(5-56) 


Then,  PQ  and  Pp  become: 


J p(y|H-|  )dy 

(5-57) 

n 

exp(-n/X-  ) 
fal 

(5-58) 

oo 

/ P(y|H0)dy 

(5-59) 
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3 exp(-n/Xr  ) (5-60) 

b0 

The  corresponding  ROC  curve  has  been  identified  as  a power  type  ROC  curve 
by  Birdsall  [25]. 


(5-61) 


(5-62) 


(5-63) 


The  ROC  curves  can  be  calculated  from  equation  (5-63)  as  long  as  the  total. 


signal-to-noise  ratio. 


, is  specified. 


A set  of  ROC  curves  with 


different  A values  are  shown  in  Figure  5.1.  The  slope  of  the  curves  in 
Figure  5.1  decreases  as  the  A value  increases.  As  shown  in  Figure  5.1 
the  slope  of  the  power  type  ROC  curves  is  always  less  than  the  slope  of 
the  ROC  curve  of  the  known-waveform  signal  case  [20].  From  equations 
(5-61)  and  (5-63),  the  parameter  A can  be  identified  as  the  total  signal- 
to-noise  ratio  plus  one.  By  using  this  relation,  we  may  look  at  Figure 
5.1  again  and  find  that  the  performance  improvement  measured  on  normal- 
normal  paper  is  not  linearly  proportional  to  the  signal-to-noise  ratio. 

In  the  high  PQ  and  low  Pp  region,  the  performance  improvement  is  less  than 
that  in  the  low  P^  and  low  Pp  region  for  the  same  signal-to-noise  ratio 
improvement.  These  observations  show  the  shortcomings  of  using  detectability 
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index  as  a performance  measure  for  the  Gaussian  signal  case. 


Example  2:  Two  Path  Signals  Detected  by  a Two-Element  Array 

In  this  example  we  assume  two  signal  paths  (M=2)  and  two  elements  (K=2) 

and  use  the  characteristic  equation  det[^  - xlj  = 0 to  solve  for  x.  Since 
2 -1 

5q  s an  1»  the  matr'x  is  equal  to  J_  . From  equation  (5-36),  we 

an 


find  the  relationship  between  \r  and  \n  as  follows: 

bl  ql 

»g,  ” 7 xq,  • ' (s-64) 

n 

Substituting  equation  (5-64)  into  equation  (5-40)  we  have: 

XG  = 1 - (5-65] 

G0  \ 

The  key  to  obtaining  \r  and  \r  is  to  solve  Xn  from  the  characteristic 

G1  R0  Q1 
equation,  det(gq  - XIJ  =0. 


(5-65) 


0,  “ *1  « <&.  + ^ v,  v,*  + as^  ^ v2*  - XI 


(5-66) 


detra, 


-X)2- 


eJuTl  ♦ o2 

5 r 


e.l^T2 


(5-68) 


The  eigenvalues  of  Q,  can  be  solved  very  easily  from  equation  (5-68). 
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2 2 2 
An  3 cT  + a.  + at  ± 

Q1  S1  s2  n 


t?  ejuTl  + ac2  eja)T2 1 


(5-69) 


The  eigenvalues  of  and  Gq  become: 


\r  = — j - 1 

Gl  an 


a?  + „ 

al  eJ"Tl  * ^ eJ“T2 

S1  s2  * 

S1  s2 

Ap  3 1 - T” 
G0  XQ 


1 


(5-70) 


(5-71) 

(5-72) 


2 2 
at  + at  ± 

a\  e-i”!  ♦ 4 eJuT2 

S1  s2 

S1  s2 

a2  + a2  + a\ 
n s-j  : 

; ± \ a2  e^UTl  + a2  eja)T2 
*2  S1  s2 

(5-73) 


The  two  non-zero  eigenvalues  of  G^  and  Gq  are  distinct  exceDt  for  the 

2 


0,  The  density  function  of  y 


rare  case  that  I a\  e^“Tl  + a eJarr2 
I s1  s2 

takes  the  form  of  equation  (5-7).  The  probability  of  detection  (Pg) 
and  the  probability  of  false  alarm  (Pp)  become: 
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o - x ~ Vx-r-  exP<-  if-  > + xr  °1,2. exD<-  if- 

Gl,l  i ,2  Gl,l  G1,2"XG1J  G1 ,2 


P„  = 


(5-74) 


pf  3 1 — ^4 exp("  r0— > + x — — exp(_  r1^ ) 

G0,l  G0,2  G0,l  G0,2  G0,l  ’0,2 


(5-75) 


where  Xr  and  Xr  are  expressed  in  equation  (5-73)  and  (5-71)  with 
u,l  bl,l 

positive  sign  and  Xr  and  Xr  with  negative  sign. 

b0,2  Q1 ,2 

In  order  to  examine  the  character  of  the  ROC  curves,  we  considered  a 
specific  case  and  plotted  the  ROC  curves  in  Figures  5.2  and  5.3.  In 
Figure  5.2  the  optimal  performance  of  a 2-element  array  for  detecting 
2-path  signals  are  plotted  for  different  total  signal-to-noi se  ratios.  The 
signal  incident  angles  are  assumed  to  be  +15°  from  the  horizontal  axis 
and  the  signal  strengths  from  both  paths  are  assumed  equal.  The  ROC 
curves  are  similar  to  the  power  type  ROC  curves  shown  in  Figure  5.1.  We 
would  also  like  to  examine  the  ROC  curves  for  2 path  signals  with  different 
signal  strengths.  The  ratios  of  signal  strength  are  assumed  to  be  1:1, 

2:1  and  10:1.  In  Figure  5.3,  three  sets  of  curves  are  drawn  on  normal-normal 
paper.  By  comparing  these  curves  we  conclude  that  the  performances  are  almost 
equal  for  low  signal-to-noise  ratios.  However,  for  high  signal -to-noise 
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ratios,  the  case  in  which  each  path  has  equal  signal  strength  gives  the 
highest  performance  among  the  three. 


Example  3 M Signal  Paths  Detected  by  a Two-element  array 

The  extension  of  the  two-path  results  in  example  2 to  the  M-path 
case  is  straight  forward.  For  the  M-path  signal  case,  0,  becomes: 


^ 3 an  — + as.|  ^1  V + + \ V 


(5-76) 


Consequently,  the  characteristic  equation  is: 


2 ......2  +02.x)2 


detfQ,  - A IJ  * (Og  +•••■*■  as 
1 M 


M 9 

i=l  1 


eJU)Ti 


(5-77) 


and  the  eigenvalues  of  Xg^  are 


M 


1 3 ac  + <*n  ± 

)i  Lu  S,  n 

1 i*l  1 


£ o\  e^i 


i^l  si 


(5-78) 


The  eigenvalues  of  Gg  and  G^  can  be  expressed  as: 

M M 


2 c2  t 

i*l  S1 

i=l  si 

(5-79) 


J 


M , M 


By  substituting  equations  (5-79)  and  (5-80)  into  equations  (5-73)  and 
(5-74),  we  may  find  PQ  and  Pp  for  the  M signal  path  case. 

Ten  signal  paths  with  equal  spatial  angular  separations  are  assumed 
to  be  within  the  range  of  +15°  from  the  horizontal  axis.  The  signal  strength 
of  all  paths  is  assumed  equal.  The  ROC  curves  of  the  2-element  array 
detector  are  plotted  in  Figure  5.4.  The  general  shape  of  the  curve  is 
similar  to  the  curves  in  Figure  5.1.  This  may  be  explained  by  the  fact 
that  there  are  only  two  non-zero  eigenvalues  and  one  of  them  is  dominant 
under  hypothesis  . Hence,  the  density  function  of  y under  is  very 
similar  to  an  expontial  curve  and  the  resulting  ROC  curves  are  similar  to 
the  power  type  ROC  curves.  The  results  shown  in  this  example  also  demon- 
strate that  it  is  not  sufficient  to  characterize  the  ROC  curves  by  using 
only  the  mean  and  variance  of  the  sufficient  statistic  y.  In  this  case, 
there  are  only  two  non-zero  eigenvalues  in  spite  of  the  number  of  signal 
paths.  The  density  function  of  y under  both  hypotheses  is  very  much  different 
from  the  Gaussian  distribution  and  cannot  be  characterized  completely  by 
only  the  mean  and  the  variance  of  y.  Besides,  one  of  the  two  non-zero 
eigenvalues  of  is  much  larger  than  another.  This  is  why  the  ROC  curves 
for  this  case,  as  shown  in  Figure  5.4,  are  very  similar  to  the  power  type 
ROC  curves. 


J 


44 
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For  any  matrix,  the  eigenvalues  can  be  evaluated  by  solvinq  the 
equation  det[Q.  - \IJ  = 0.  However,  if  the  dimension  of  the  matrix  is 
larger  than  4,  the  eigenvalue  cannot  be  solved  for  analytically.  Even  for 
a (3x3)  matrix  case,  obtaining  the  eigenvalues  can  be  quite  tricky. 

Example  4 will  show  a case  with  two  signal  paths  and  a 3-element  array. 

Example  4 Two  Signal  Paths  Detected  by  a Three-element  Array 

We  assume  M=2  and  K=3  in  this  example,  where  M is  the  number  of  paths 
and  K is  the  number  of  array  elements.  The  covariance  matrix  becomes: 

3,  ■ ^ I * V * a\2  % % (5-81) 

Then, 

9q  ' =^CTn  “ X)-  + as1  -1  -1  + as2  ^2  ^2 


an  ‘ X 

o 

o 

0 

1 eJWTl  e2jarrl 

o 

^-x 

o 

e_ju,Tl  ! <>Ti 

o 

o 

an  ~ X 

e-ja)T1ae-ja)T1  i 
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+ a. 


eJwt2  e^^'UT2 


rju)T2 


ej“T2 


■2jajT2  a-JWT, 


e JW12 


(5-83) 


and. 


detfa-i  ‘ x-^ 


<4,  * °l2  * 4 - *>*  * <N>2  *2  * A1  A2* 


-<«?,  * ^ * 4 - *><V2*  * 2AiAi*> 


(5-84) 


where 


A = a?  eJUTl  + a?  ejaJT2 
1 Sl  s2 


(5-85) 


and 


A.  = ct?  e2'iwTl  + a2  e2jwT2 


2 us 


(5-86) 


From  equation  (4-29)  we  know  that  the  matrix 


■1 


9<>  - ^ 3Ct  v-,  v,  + c2  VgVg  - (yg  v^jvgv,  + (v,  v^Vg  ] (5.87) 


has  rank  2.  We  also  know  that  both  Gq  = (Q^1  - Qq"1)  and 


G-j  = (Pq*1  - Q.q"1)  have  rank  2.  This  implies  that  one  of  the  solutions 
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to  the  characteristic  equation  of  G^  and  Gj  has  to  be  zero.  In  other 
words,  one  of  the  eigenvalues  of  should  satisfy  the  following  two 
equations. 

c2 
'G 


X.  3 1 - r11-  3 0 


Jo 


^i 


-1=0 


(5-88) 

(5-89) 


We  found  that  Xg  should  be  equal  to  a2  to  satisfy  both  equation  (5-38) 
and  equation  (5-89).  This  eigenvalue,  Xg  = a2  , will  be  verified  with 
the  following  equations. 


det[^  - XQ  Ij  3 detta,  * an  U 


= detfa^  v,v,*  + °l2  VgVg*] 


(5-90) 

(5-91) 


Since  the  rank  of  a2  v^v-,*  + a2  VgVg*  is  two  and  the  dimension  of  the 

2 

matrix  is  three,  we  can  easily  see  that  det[^  ~ crn  U = 0,  This 
shows  that  Xn  3 a2  is  indeed  a solution  to  the  characteristic  equation 
of  a,.  The  next  step  is  to  reduce  the  3rd  order  equation  by  dividing 
detfQq  “ XIJ  by  (a2  - X).  After  some  manipulations  we  obtained  a simple 
second  degree  equation  in  x. 


il 


47 


det[Q_^  - Xlj  2 2 2 2 2 9 9 9 

O * 3(af  + of  )Z  + 3(o2  + of  )(o2  - X)  + (of  - X)2 


°n  • X 


S1  s2 


S1  s2  n 


"[^2^2  + 2^1^1  3 s 0 


(5-92) 


where  and  A2  are  defined  in  equations  (5-85)  and  (5-86). 

We  may  solve  the  second  degree  equation  to  get  xn  . 

gl 


<°n  * X> 


■3(aj  + o2  ) 
S1  S2 


1 J '3(^1  * ^2>2  * 4(^*  * 


(5-93) 


so. 


3(os  + °s  ) 
X - 2 + 1 2 

\ ' an  2 


-i>J  "3(os1 


+ a^)2  + 4(A2A2*  + 2A1A]*) 


(5-94) 


"here  ’3  (°S,  + »s2)2  + «(V2*  * 2Vl‘)  • 9 ”s,  °<2  - 6 ^ ,*[ 

^ - j2a)Ti  + i2(i)T2  -j2uT2  + j2(»)t^  ^ - jtjJT  1 + jtOT  2 ju)t^-j(UT2^J 
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Substituting  equation  (5-94)  into  equations  (5-70)  and  (5-72),  we  can 

get  \r  and  Xr  . 

G0  G1 

To  show  the  detection  performance  of  the  3-element  array,  we  use 
the  same  two  signal  paths  as  in  example  3.  The  optimal  ROC  curves  for 
several  signal-to-noise  ratios  are  shown  in  Figure  5.5  for  a 3-element 
array.  We  notice  that  there  are  only  2 non-zero  eigenvalues  because 
there  are  only  two  paths.  In  general,  the  number  of  non-zero  eigen- 
values is  no  larger  than  min(K,M),  where  K is  the  number  of  array  elements 
and  M is  the  number  of  signal  paths.  In  Figure  5.5,  the  shape  of  the  ROC 
curves  for  a 3-element  array  is  similar  to  that  in  Figure  5.2,  and  is 
apparently  not  like  normal  ROC  curves  which  would  be  the  result  for  the 
case  with  Gaussian  distributed  sufficient  statistics  under  both  hypo- 
thesis. By  using  the  equation  A = 1 + K-SNR,  we  also  recognize  that  the 
slope  of  the  ROC  curves  in  this  example  is  slightly  higher  than  that  for 
the  power  type  ROC  curves  shown  in  Figure  5.1. 

In  order  to  see  the  performance  improvement  by  adding  hydrophones 
to  the  array  receiver,  we  make  the  performances  comparison  among  arrays 
with  different  numbers  of  elements  in  the  next  section. 

5.3  Performance  Comparison  between  Optimal  Array  Processors  with 

Different  Numbers  of  Elements 

In  the  previous  section,  we  have  derived  performance  evaluation 
methods  analytically  for  some  special  cases.  However,  we  will  not  solve 
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every  problem  by  using  the  analytical  method  because  it  is  too  tedious 
and  error  prone.  A numerical  method  for  finding  the  eigenvalues  of  a 
complex  matrix  is  shown  in  Appendix  A.  By  employing  the  software  sub- 
routine package  EISPACK,  a computer  program  for  evaluating  the  detection 
performance  numerically  has  been  developed  to  find  the  performance  of 
arrays  with  more  than  3 elements.  To  make  sure  of  the  correctness  of 
the  computer  program,  we  have  checked  the  numerical  results  with  the  ROC 


curves  obtained  from  the  analytical  method  for  several  cases. 

We  use  the  case  with  two  signal  paths  coming  in  at  +15°  from  the 
horizontal  axis  as  an  example  to  show  the  detection  performance  improve- 
ment due  to  adding  hydrophones  to  the  array.  In  Figure  5.6  and  Figure 
5.7,  performances  of  arrays  with  one,  two,  and  three  elements  are  plotted 


for  different  values  of  signal-to-noise  ratio  per  element,  SNR  = 


2.2 

cs  °s 
S1  s2 


2 2 

where  a and  a.  are  signal  power  of  path  1 and  path  2 respectively. 
S1  s2 


The  ROC  curve  has  a slope  less  than  tan  45°  on  the  normal -normal  paper.  In 
general,  the  slope  of  the  ROC  curves  is  smaller  for  high  signal-to-noise 
ratios  than  for  low  signal-to-noise  ratios.  We  also  notice  that  the 
performance  improvement  obtained  by  adding  one  hydrophone  to  a 2-element 
array  is  less  than  that  obtained  by  adding  one  hydrophone  to  a single-element 
array  for  SNR=10.  However,  for  low  signal-to-noise  ratios  such  as  SNR  of 
1 and  2,  improvements  from  adding  one  hydrophone  to  a 2-element  or  a 3-element 
array  are  nearly  the  same.  This  phenomenon  can  also  be  seen  by  comparing 
the  solid  lines  in  Figures  5.8  and  5.9.  In  Figure  5.8,  the  ROC  curves  for  SNR=1 
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and  K=1 , 2,  3 and  4 are  parallel  with  equal  spacing.  On  the  other  hand, 
the  ROC  curves  in  Figure  5.9  for  SNR=10  and  K=1 , 2,  3 and  4 have  unequal 
spacing.  Apparently,  the  performance  improvement  due  to  adding  one  hydro- 
phone to  an  array  depends  on  the  array  factor  as  well  as  the  signal -to- 
noise  ratio  per  element.  In  general,  the  performance  improvement  is  less 
for  high  SNR  than  that  for  low  SNR.  This  shows  a difference  between  the 
performance  of  detecting  multipath  Gaussian  signals  and  that  of  detecting 
multipath  known  waveform  signals  [28]  in  which  the  performance  improvement 
due  to  adding  hydrophones  is  independent  of  signal-to-noise  ratio. 

It  is  also  interesting  to  compare  the  performance  characteristics 
of  array  detectors  with  that  of  a single  phone  detector.  We  look  into  the 
question  of  how  much  performance  improvement  an  array  can  gain  over  the 
signal  phone.  To  answer  this  we  compare  the  performance  of  a single  hydro- 
phone detector  with  the  performance  of  the  array  detector  with  K elements. 

In  Figure  5.8  we  compare  the  performance  of  a single  hydrophone  detector 
with  K times  the  signal-to-noise  ratio  with  the  performance  of  K-element 
array  detectors  for  K=1 , 2,  3,  4 and  10,  and  SNR=1 . We  observe  that  the 
performances  of  these  two  cases  are  nearly  the  same  for  small  (K*SNR)  values. 
But,  for  large  (K«SNR)  values  the  performance  of  an  array  detector  is  better 
than  that  of  the  single  hydrophone  case.  In  Figure  5.9  we  make  a similar 
comparison  for  SNR»10.  The  curves  in  Figure  5.9  show  dramatic  performance 
improvement  for  the  array  detectors.  By  referring  to  the  structure  of  the 
optimum  array  detector  in  which  signals  from  each  path  are  spatially  separated 
from  other  paths  by  beamforming  action,  we  may  conclude  that  the  array 
detectors  not  only  receive  a total  output  signal-to-noise  ratio  K times 
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higher  than  that  of  the  single  phone,  but  also  provide  the  capability  of 
separating  multipath  signals  spatially.  The  spatial  discriminating 
capability  of  the  array  detector  does  improve  the  performance  considerably 
for  high  signal-to-noise  ratio  cases. 

By  inspecting  the  slopes  of  the  curves  in  Figure  5.8  and  5.9,  we 
notice  the  difference  in  shape  between  the  array  and  single  hydrophone 
ROC  curves.  This  can  be  explained  by  looking  at  the  density  function  of 
y.  As  discussed  before,  we  know  that  for  a single  hydrophone  the  density 
function  of  y is  exponentially  distributed.  On  the  other  hand,  the  density 
function  of  y for  the  array  detector  depends  on  two  non-zero  eigenvalues 
for  various  K values  and  signal-to-noise  ratio  as  is  shown  in  Table  I. 

From  the  table  we  observe  that  in  some  cases  two  eigenvalues  are  almost 
identical  and  in  other  cases  one  eigenvalue  is  larger  than  the  other.  For 
the  K=2  cases,  there  exists  a dominant  eigenvalue  which  is  due  to  the  fact 
that  the  size  of  the  array  is  too  small  to  discriminate  two  signal  paths. 
Hence,  the  resulting  ROC  curves  have  shapes  similar  to  that  of  the  power 
type  ROC  curves,  which  is  the  character  of  the  ROC  curves  for  a single 
hydrophone  detector.  For  the  K=4  cases,  the  width  of  the  mainlobe  of  a 
4-element  array  is  about  60°  which  is  about  two  times  the  angular  separation 
between  the  two  signal  paths.  In  other  words,  the  4-element  array  may 
point  a beam  to  one  path  and  null  out  the  other  path  simultaneously.  As 
a result,  the  eigenvalues  for  this  case  are  almost  identical.  Consequently, 

p 

the  distribution  of  y is  the  x distribution  with  4 degrees  of  freedom. 

The  noticeable  difference  on  ROC  curves  between  the  exponential  distribution 
2 

and  the  x distribution  with  4 degrees  of  freedom  is  the  slope  of  the  ROC 


A 


52 


curves.  The  curves  of  the  former  case  have  smaller  slopes  than  that  of 
the  latter  case. 

By  investigating  the  beam  pattern  of  a K-element  array  [27],  we  found 

that  the  main  beam  width  is  approximately  equal  to  1 80°/ (K-l ) . For  K=2 

and  K=3  cases,  both  signal  paths  are  in  the  mainlobe  of  the  beam  pattern, 

and  for  K=4  and  K=10  cases,  one  path  is  in  the  mainlobe  and  the  other  is 

in  the  si  delobe.  Consequently,  the  eigenvalues  of  for  the  K=4  and 

K=10  cases  are  close  to  each  other  while  the  \r  's  for  K=2  and  K=3  are 

G1 


far  apart.  The  dependency  of  the  eigenvalues  on  the  array  beam  pattern 
will  be  discussed  further  in  the  next  section  in  which  we  will  investigate 
the  performance  of  suboptimal  processors  such  as  single  and  double  beam- 
formers  . 


L. 


A 


TABLE  I 
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5.4  Detection  Performance  of  Suboptimal  Array  Beamformers 
A.  Single  Beamformer 

A single  beamformer  which  points  to  one  of  the  M signal  paths  and 
a double  beamformer  which  points  to  two  signal  paths  are  investigated 
in  this  section.  The  trade-off  between  the  complexity  of  the  processor 
structure  and  the  detection  performance  will  be  evaluated. 

A single  beamformer  is  a suboptimal  array  which  ignores  the  existence 
of  the  multipath  structure  and  points  a beam  to  the  direct  path.  The 
simplicity  of  the  processor  structure  makes  it  popular  in  practical 
applications.  The  sufficient  statistic  y of  the  single  beamformer  may  be 
expressed  by 


★ ★ 

y = R u u R 


(5-95) 


where  £ is  the  pointing  vector  of  the  single  beamformer. 
Then,  by  definition. 
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and 
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where 


(5-100) 
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It  is  quite  easy  to  find  the  eigenvalue  of  Gq  from  the  next  equation. 


Gq  - * n - (-  -) 


3 Ka  u 
n — 


(5-101) 

(5-102) 


To  find  the  eigenvalues  of  is  not  obvious.  We  should  use  the  matrix 
theorem  that  "the  eigenvalue  of  G^  is  equal  to  that  of  G^".  By  observa 
tion,  we  found  the  eigenvector  and  eigenvalue  of  from  the  following 
equations. 
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*T 


i aVa! *T>  * 2 
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(5-103) 
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Hence, 
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(5-105) 


The  ROC  curve  for  this  case  is  apparently  a power  type  ROC  curve  and  can 
be  expressed  by 


(PD> 


A 


(5-106) 
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where 


(5-107) 


n 9 
i ♦ 


^ U| 


to: 


(5-108) 


The  second  term  in  equation  (5-108)  is  the  summation  of  signal  energy 
weighted  by  the  beam  pattern. 

In  order  to  compare  the  single  beamformer  performance  with  the  optimal 

processor  performance,  two  signal  paths  are  chosen  coming  at  +15°  from  the 

horizontal  axis.  The  values  of  A for  K = 2,  3,  4 and  10  are  2.4723,  2.6313, 

3.0029,  6.1121  respectively  in  the  cases.  We  notice  that 

S1  s2  n 

the  value  of  A can  be  approximated  by  A = 1 + — when  K>4.  This  is  due 


to  the  fact  that  the  beamwidth  of  the  main lobe  for  K>4  is  less  than  60°. 
Hence,  the  second  signal  is  in  the  sidelobe  of  the  beam  pattern  and  does 
not  contribute  too  much  signal  energy.  In  other  words,  the  single  beam- 
former  array  nulls  out  the  second  signal  source  when  the  size  of  the  array 
is  large.  Consequently,  it  has  the  same  performance  as  that  of  the  single 

K cr? 


2 * 


path  case  in  which  the  ROC  curves  are  characterized  by  A = 1 + 
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The  ROC  curves  of  the  single  beamformer  are  plotted  in  Figure  5.10 
and  Figure  5.11  for  different  signal-to-noise  ratios.  By  comparing 
these  ROC  curves  with  that  of  the  optimal  detectors,  we  conclude  that  the 
performance  loss  due  to  not  taking  account  of  the  second  path  increases 
as  the  number  of  array  elements  increases.  This  indicates  that  the  larger 
the  array  size,  the  more  important  the  knowledge  about  the  spatial  structure 
of  the  si gnal . 


B.  Double  Beamformer 

If  we  want  to  retain  the  beamforming  structure  and  improve  the  per- 
formance, then  the  double  beamformer  is  our  next  choice  for  detecting 
two  path  signals.  The  sufficient  statistic  y is  described  in  equation 
(5-109). 
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2*2  * 

eigenvalues,  X0B,  of  matrix  (as  VqV^  +a  s v^Vg  ),  because  XG  and 
XG  may  be  expressed  In  terms  of  XBB  as  shown  in  the  following  equations. 


XDB 


(5-114) 


XG1  = an  X0B  + XDB  (5-115) 

Eigenvalue  XDB  may  be  evaluated  by  both  analytical  and  numerical  methods. 

A simple  two-element  case  will  be  shown  as  an  example  ^or  the  analytical 
method. 
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Substituting  equation  (5-118)  into  equations  (5-114)  and  (5-115),  we  may 
find  Aq  and  \g  . The  rest  of  the  performance  evaluation  procedure  is 


similar  to  that  of  the  examples  in  the  previous  sections. 

The  structure  of  the  double  beamformer  can  be  realized  as  two  beam- 
formers  followed  by  envelope  detectors  as  shown  in  Figure  5.12.  Comparing 
Figure  5.12  with  Figure  4.1,  the  only  difference  is  that  the  cross  correlated 
terms  are  missing  in  Figure  5.12.  The  performance  degradation  due  to  the 
lack  of  the  correlated  terms  is  the  main  topic  in  this  subsection. 

The  performance  of  the  double  beamformer  has  been  calculated  by  using 
numerical  programs  for  two,  three,  four,  and  ten  element  arrays.  The  ROC 

curves  are  shown  as  the  dashed  lines  in  Figure  5.13  and  Figure  5.14.  On 
the  same  Figures  the  ROC  curves  of  the  optimal  detectors  are  also  shown. 

It  is  interesting  to  see  that  the  performance  of  the  double  beamformer  is 
very  close  to  optimum  for  K>4.  This  can  be  explained  by  the  fact  that  one 
of  the  two  signal  paths  is  in  the  sidelobe  of  the  beam  pattern  and  consequently 
the  cross  correlated  terms  in  the  optimal  array  are  small  compared  to  the 
beamformer  output  power.  Hence,  the  cross  correlated  terms  may  be  ignored 
when  the  array  size  is  large  enough  such  that  the  main  beam  width  is  less 
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than  two  times  the  angular  separation  of  the  two  signal  paths.  In 
Figure  5.13  we  observe  that  the  performance  of  the  double  beamformer 
is  very  close  to  optimum  for  K = 2,  3,  4 and  10  at  a signal -to-noise 
ratio  of  1.  But,  in  Figure  5.14  the  performance  degradation  is  severe 
at  K = 2 at  a signal -to-noise  ratio  of  10.  This  indicates  that  the  cross 
correlated  term  plays  an  important  role  only  for  cases  with  high  signal- 
to-noise  ratio  and  small  array  sizes. 


CHAPTER  VI  SUMMARY 


This  report  has  proposed  vertical  arrays  for  estimating  parameters 
of  multipath  signals  and  for  detecting  multipath  Gaussian  signals.  The 
spatial  discriminating  capability  of  vertical  array  processors  can  be 
used  to  isolate  individual  signal  paths  in  an  underwater  acoustic  channel. 
The  amplitude  and  phase  estimates  of  signals  from  each  individual  path 
are  obtained  and  compared  to  that  for  a single  hydrophone  case.  The 
signal-to-noise  ratio  of  each  individual  path  is  shown  to  be  independent 
of  the  relative  signal  phases  of  other  paths.  This  suggests  that  the 
vertical  array  processor  may  be  useful  in  the  deep  fading  problem  for 
minimizing  multipath  signal  cancellation.  Two  approaches  are  taken  to 
estimate  the  signal  parameters.  In  Chapter  III,  the  signal  parameters 
are  assumed  non-random  and  the  maximum  likelihood  estimates  (MLE)  and 
their  performances  are  shown.  It  is  clear  that  if  the  number  of  array 
elements  is  less  than  the  number  of  signal  paths,  then  the  maximum  likeli- 
hood ratio  estimate  is  not  unique.  In  Chapter  IV,  the  signal  parameters 
are  assumed  to  be  Gaussian  random  variables  with  covariance  matrix 
and  the  mc-imum  a posteriori  probability  estimates  (MAP)  of  multipath 
signal  parameters  are  employed.  With  the  a priori  knowledge  about  the 
signal  in  terms  of  the  covariance  matrix  the  MAP  estimate  has  a unique 
value  even  for  the  case  where  the  number  of  array  elements  is  less  than 
the  number  of  signal  paths.  The  linear  relationship  between  MAP  and  MLE 
estimates  is  presented  to  show  the  incorporation  of  the  a priori  knowledge 
of  the  signal  parameters. 


The  second  part  of  this  report  is  devoted  to  the  detection  of 
multipath  Gaussian  signals  in  Gaussian  noise.  Like  most  Gaussian  signal 
detectors  the  optimal  array  detector  has  a quadratic  structure.  A 
two-path  signal  example  shows  that  the  optimal  array  not  only  points  two 
beams  in  both  signal  directions  but  also  correlates  both  beams  to  form  a 
cross  correlation  term,  even  when  multipath  random  signals  are  assumed 
independent  from  path  to  path.  In  Chapter  V,  the  significance  of  the 
two-beamformer  structure  and  the  cross  correlation  term  are  investigated 
in  terms  of  detection  performance.  An  analytical  method  called  the 
"Eigenvalue  method"  is  developed  for  evaluating  the  performance  of  any 
array  signal  detector  which  has  a quadratic  sufficient  statistic.  The 
eigenvalue  method  starts  with  the  characteristic  function  of  the  sufficient 
statistic  under  both  hypotheses  and  then  transforms  the  characteristic 
function  to  the  probability  density  function.  The  probability  of  false 
alarm  and  probability  of  detection  may  be  calculated  analytically  once 
the  eigenvalues  of  the  covariance  matrices  of  the  observed  data  under  both 
hypotheses  are  calculated.  Although  we  may  find  the  eigenvalues  of  any 
matrix  with  dimension  less  than  four  by  an  analytical  method,  the  calcula- 
tion of  eigenvalues  is  performed  mostly  by  a numerical  method  in  Chapter  IV. 
For  matrices  with  dimensionality  less  than  three,  both  methods  are  demon- 
strated. 

Performance  comparisons  were  made  between  optimal  arrays  with  different 
number  of  elements.  By  use  of  the  two-path  signal  case  as  an  example,  we 
found  that  the  performance  improvement  obtained  by  adding  one  hydrophone 
to  a 2-element  array  is  less  than  that  of  adding  one  hydrophone  to  a 
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single-element  array  for  high  signal-to-noise  ratios  around  SNR  = 10. 

But,  for  low  signal-to-noise  ratios  such  as  one,  two  and  four  the  performance 
improvement  for  both  cases  are  nearly  the  same.  This  shows  the  difference 
between  the  performance  of  detecting  Gaussian  signals  and  the  performance 
of  detecting  known  waveform  signals.  It  is  well-known  that  in  the  latter 
case  the  detectability  index  is  the  product  of  two  functions,  one  depends 
on  the  signal-to-noise  ratio,  another  depends  on  the  array  factor.  In 
other  words,  the  performance  improvement  in  terms  of  detectability  index 
by  adding  one  hydrophone  to  an  array  is  independent  of  the  signal-to-noise 
ratio.  The  phenomena  of  the  former  case  may  be  explained  by  the  property 
of  power  type  ROC  curves  which  are  the  type  ROC  curves  that  result  in  the 
Gaussian  signal  detection  problem.  The  family  of  power  type  ROC  curves 
shows  that  the  performance  improvement  due  to  the  increase  in  signal-to- 
noise  ratio  starts  to  decrease  after  the  parameter  A becomes  larger  than 
about  21.  During  the  process  of  comparing  the  ROC  curves  for  the  Gaussian 
signal  array  detectors  with  the  power  type  ROC  curves,  we  observed  the 
difference  in  the  slopes  of  ROC  curves.  This  led  to  a comparison  between 
the  performance  of  a k-element  array  with  that  of  a single  phone  detector 
with  k times  the  signal-to-noise  ratio.  For  low  total  output  signal-to- 
noise  ratios  in  the  range  of  one  to  four  the  difference  in  performance 
between  these  two  cases  is  insignificant.  However,  for  high  signal-to- 
noise  ratios,  (i.e.  (k-SNR)>10) , the  array  detector  is  better  than  the 
single  phone  detector,  for  the  same  total  signal-to-noise  ratio.  This 
may  be  attributed  to  the  spatial  discriminating  capability  of  the  array. 


The  trade-off  between  the  processor  complexity  and  performance  is 
also  shown  in  Chapter  V by  comparing  the  performance  of  an  optimal  array 
detector  with  that  of  conventional  beamformers.  Again,  the  two-path  case 
is  used  to  demonstrate  the  performance  difference  between  optimal,  single 
beamformer  and  double  beamformer  arrays.  A single  beamformer  array  which 
points  a beam  in  one  of  the  two  signal  directions  performs  poorly  com- 
pared to  the  optimal  array.  Although  the  structure  of  the  single  beam- 
former  is  only  one  branch  of  the  optimal  processor,  the  performance  of 
a 10-element  single  beamformer  array  is  equivalent  to  that  of  the  optimal 
array  with  only  4 elements  at  SNR  = 1. 

Another  suboptimal  array  processor  which  points  two  beams  to  two 
signal  paths  and  sums  up  the  envelope  detector  outputs  is  also  investigated. 
The  performance  comparison  shows  that  if  the  aperture  of  the  array  is  large 
enough  such  that  the  second  signal  direction  is  not  in  the  mainlobe  or  any 
significant  si  delobe  of  the  array  pattern,  the  difference  between  the 
suboptimal  array  and  optimal  array  in  performance  is  very  small  for  all 
signal-to-noise  ratios.  For  high  signal-to-noise  ratios  and  small  arrays 
the  performance  difference  between  the  suboptimal  and  optimal  arrays  is 
noticeable.  This  finding  implies  that  the  cross  correlation  term  in  the 
optimal  processor  is  not  very  significant  if  the  array  aperture  is  large 
enough  to  isolate  signal  paths  or  the  signal-to-noise  ratio  is  low. 

The  eigenvalue  method,  introduced  in  this  report,  is  a powerful  tool 
for  evaluating  performance  of  array  detectors.  Since  the  detection  per- 
formance and  the  processor  complexity  are  two  major  factors  to  consider 
in  the  implementation  of  array  detectors,  the  trade-off  of  these  factors. 
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as  shown  in  this  report,  gives  insight  into  the  problem  of  detecting 
multipath  Gaussian  signals  in  Gaussian  noise. 
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APPENDIX  A 

COMPUTER  PROGRAMS 

The  computer  programs  for  evaluating  eigenvalues  of  matrices 
Gq  and  G^ , and  for  calculating  Pg  and  Pp  are  shown  in  this  appendix. 

The  subroutine  calling  sequences  are  shown  in  the  block  diagrams.  There 
are  two  types  of  methods  for  evaluating  the  eigenvalues  of  a complex 
matrix.  One  is  the  analytical  method  which  is  limited  by  the  dimen- 
sionality of  the  matrix.  It  is  not  trivial  to  find  the  performance  of 
array  detectors  with  more  than  three  elements  by  the  analytical  method. 
Another  method  is  a numerical  one  which  can  be  used  for  arrays  with 
any  number  of  elements. 

(1 ) Analytical  Method 

A typical  calling  sequence  for  the  analytical  method  is  shown  below. 


Subroutine  "EIGEN"  calculates  the  eigenvalues  of  matrices  Gg  and  G^ , and 
subroutine  "PDPF"  calculates  the  probability  of  detection  (PQ)  and  the 
probability  of  false  alarm  (Pp).  The  analytical  algorithms  used  in 


A-l 


subroutine  "EIGEN"  depend  upon  the  number  of  array  elements,  the  number 
of  signal  paths,  and  the  sufficient  statistic.  The  mathematical  descrip- 
tions of  these  algorithms  are  described  in  the  main  text.  In  the  follow- 
ing computer  program  listings,  a complete  program  for  optimal  array 
detectors  for  K = 2 and  M = 2,  several  "EIGEN"  subroutines,  and  a complete 
program  for  optimal  array  detectors  for  M = 10  and  K = 2 are  included. 
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C THIS  IS  ThE  ma[h  kGh  FOR  CALCULATING  hope  a r ANALIUCAL  "ETriOu. 

C T-tlS  PGM  «ILL  CALL  SUoRUUT  I nE  E I Gt  N whILh  *ILL  CALCuAIE  E iGt  N v ALUfcS 
C OF  G0,Gl.  ThE  M a I i'i  ALSU  CALL  aUbRuUri'xE  PlfPF  TO  GET  Pu  AN  u pF 
C w£  aSSSume  Us.5  WAVELENGTH 

DIMENSION  EL  l A ( J u ) ,Pu( 1 0 ) #PF ( 10  ) 

WHITE  (o.lOn 
REA0(5*l0u)  APHI,aPhIP 
PHlsAPHI*i.l4ib/lbO. 

PhIPsAPHlP*i.l416/lbO. 
wHI  TE  Cb#  1 o2 ) 

HE AO ( b f 100)  VARSl# VARS2#VaHn 
wRI  TE(b, 103) 

REAO(b»  100 ) (ELTA ( I ) , I = i , 10  ) 

100  F0RMAT(F14.0) 

101  FuRMAT ( 1X,24HTYPE  In  APh I , aphj P , F 1 4 .0  ) 

101  FjRfiA i ( ix,2ohttpe  in  io  threshold  value) 

102  Format ( IX, 30HTYPE  in  v akS  1 r v arS2, v ahn, F i u . u ) 

CALL  EiGENiPHl,PhIP,VARSl#vARS2*vAKN,Gil,Gl2,G0l.G02) 

CALL  POPF(ElTA,G1 l#G12#G0l»G02#HO,PF) 
wRITE(6,2uO)  Gll,G12,GOl,G02 
WRITE (6, 300)  CPU(J),PF(J),J=1, 10) 

300  FORMAT (2X,2F12.o) 

200  F0RMATC2X#2Flv).4) 

STOP 

E^D 

SUdRUUT InE  P0PF(ELTA,G1 1,G12,G01,G02,P0#PF) 

DIMENSION  ElTa(10),PD(10)»PF(10) 

DU  10  K=l,iO 

PD(K)=G1  1*ExPC-ELTA(K)/G1  1 )/(GU-G12)+Gl2*ExPC-ELTA(K)/ 
l G12)/(Gl2-Gll) 

PF (K)=GOl *ExP(-£LTAtK) /G01 ) / ( GO  1 -b02 ) *Gi)2 *E *P t -EL T A t* ) / 

1 G02)/IG02-G0t ) 

10  CONTINUE 

RETURN 

End 

SUBROUTINE  EIGEN(PHI,PHIP,  VaRSI*  VARS2#  VaRn,GU.G12#G01,G02) 

C THIS  IS  FUR  Ks2,M=2  OPTIMAL  CASE  (AnAlITICAL  -iEThuD). 

C wE  ASSSUME  o* .5  WAVELENGTH 

CCl  = 3,  t4lo*COS(PHl ) 

CC2=3. 1416*C0SCPHIP) 

COMPLEX  ai,a2#C1,C2 

ci=cmplxco.,cci) 

C2=CMPLXlO.#CC2) 

A1=CEXP(C1)  . 

A2=CEaP(C2) 

R=CAHS( VARSI*A1 >VARS2*A2) 

Gi  l = (VARSl fVARS2+R)/VARN 
Gl2s(VARSltVARS2-R)/VARN 

G01=(VARSltVARS2+R)/CVAR31>VARS2tVARN+R) 

G02=( VARS1+VARS2-H)/ ( VARSl^v ArS2+vAkn-h ) 

RE  TURN 
END 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
RiOltf  COf  YEURNISHED  XODDC  


C This  IS  SuBHUuT  [lit  EIGEn  PGR  QOuBlE  BtAHp'uRMER 
C USE  This  *IfH  main  PGm  AMU  SUBROUTINE  PuPP 

SUHROu  T I N£  EIGEN  (PH  I,  PH  IP,  VAWSl»  VaRS2,  VAKN,G1 l,bl2,G0l,G02) 
c This  IS  pur  SuBuPTIhal  BE  AMP  OWmER  P£RPJnManlE  K=2,ms£ 

CUMPlEx  C1,C2 

CI=CMPtX(0.,i, 1416»CUSCPHl) ) 

C2=CMPLxCo.,3.l4lo*COS(PHlP) ) 

OA1SVAHS1 tVARS2fCABS( VARS1 *LE *P l C 1 ) ♦ V ARS2 *C E XP ( C2 ) ) 
i)A2  = 2,*l  vARSl»VARS2)-0Ai 
GO l sV AHn*OA 1 
G02sVARN*QA4 
G1  1svAhN*DAUOA1*OA1 
Gl2bVARN*UA2+DA2*DA2 

return 

END 


SUBROUTINE  EIGEN (PHI, PH  IP, VaRSI. VARS2, VARN.Gl l,Gl2,G0l  ,G02) 

C This  IS  The  SUBROUTINE  EIGEN  FOR  CALCULATING  EIGENVALUES  FOR  Ms2,K5i 
C This  IS  FOR  Ks3,M=2  OPTIMAL  CASE  (ANALITICAL  METHOD). 

C USE  This  wITH  main  PGM  AND  subroutine  pdpf 
COMPLEX  CC I , CC2 

CC2sCMPLXCO.,6.2832*(COS(PHI)-CUS(PHlP) )) 
CClsCMPLX(0.»3.l41b*(CQS(PHl)-COS(PHlP) ) ) 

As9,*(VARSl »*2+VARS2**2)-b.*VARSl *VARS2*4, *v AnSl *vARS2*12. 

1 *REAL(CEXP(CC2))P4.*REAL(CEXP(CC1))) 

Q] lsvARNt3.*CVARSltVARS2)/2.tSQRT(A)/2. 

Q1 2 = V4RNi-3.*C VARSI +V ARS2)/2.-SGRf  (A)/2. 

GllsQU/VARN-1. 

G12=012/VARN-1. 

G01s1,-VARn/Q11 

G02*U-VARn/UI2 

RETURN 

END 


C This  PRG  IS  FOR  MstO  AMO  K = 2 OPTIMAL  case 
C This  IS  THE  PGM  To  CALCULATE  PD  And  pf 

DIMENSION  ELTA(1U),PQ(10),PF(  10),PHIUU) 

M=  1 0 

DO  5 KK=1,5 

PHI  (KK)s(90.-FLQAr(*K)  *3 . } *3 . 1 4 1 6/ 1 oO  . 

KKlstO-KK+l 

PHl(KKl)s(90.+FLUAT (KK) *3.)*i, 1416/1B0. 

S CONTINUE 

wWITECb, 102) 

HE AD ( B » 100)  VAR'N 

102  FORMAT(1X,1‘2hTYPE  IN  VARN)  THIS  PAGE  IS  BEST  QUALITY  PRAC1ICABL1 

wRITE(6,lo3)  Fim  G0|PI  E(URiilSHED  to  ldq 

REAOCS, 100)  CELTA(I),Isi,iO) 

100  FQRM A T ( F 1 4 . 0 ) 

VAHS=l ,/FLUA  T(M) 

103  FORMAT! IX, 2oHTYPE  In  10  THRESHOLD  VALUE) 

CALL  EIGtN(M,PHl,VARS,vARN,Gll,Gl2,G0l,G02) 

CALL  P0PF(ELTA,Gtl,G12,G01,G02,PD,PF  ) 
wRl  TE (b, 200  ) G 1 1 , G 1 2 , GO  1 , G0  2 
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»Rl TE 300  ) lPo( J ) / PF  C J ) , J=1 , lu  ) 
iOu  FjRmA  [ (2*#2|- l2.o) 

200  forma  r ( 2X  , 2F  1 U .4  ) 

s top 

EnO 

SUoRijUT  I it  PoPF  (EL  T A,Gl 1 ,G12,GU1 ,Gu2, Pu#PF  J 
D I mens  I uN  EL  T A ( l 0 J , Pi)  (l  u ) , PF  ( 1 0 ) 

Du  10  Ksl, 10 

IPDCisJsGi  l*fcxPl-tLTA(KJ/Gl  l)/(Gll-Gl2)tGl2*E*P(-tLfA(Kj/ 

1 G 1 2 ) / ( G 1 2-G 1 t ) 

PFU)=Gul *EaP(-ELT  A(k)/GOi  )/ ( GO  1 -U02 ) ♦U02*EaP ( -fe L T A ( K i / 
l G02)/(G 02-G01) 

10  CONTINUE 

rETuRn 
EnO 

subroutine  eige.n(r,phi, vars#var.n,gi i#i»12#gui»G02) 

DIMENSION  PnIU0),CU0),A(lQ) 

COMPLEX  C * A , 8 
BsCmPLX ( 0 . * 0 • ) 

I DO  10  1=1, M 

C(I)=CMPlX(0#,  3.Ulb*COSCPHICI))) 

A ( i ) =CEXP (C ( I ) ) 
a=d*A(I)*VARS 
10  CONTINUE 

RsCaBS(o) 

AMsFLJA  T(M) 

Gii  = (am*vars-*-«)/varn 

G12=(AM*wa«S-R)/VARN 

G01  = (A!i*VARStP)/lVARN^AH*VARS  + K) 

Go2=(AM»VARS-H)/CVARN*AM*VAHS-«) 

RETURN 

END 

' 
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(2)  Numerical  Method 


A subroutine  package  called  "EISPACK"  [24]  is  employed  in  the 
numerical  method.  The  calling  sequence  is  described  in  the  block  diagrams 
below  for  both  the  optimal  case  and  the  double-beamformer  case. 


Figure  A.  1 Block  Diagram  for  Optimal  Case. 

r 


Figure  A.  2 Block  Daigram  for  Double  Beamformer  Case 

The  following  program  listing  includes  the  main  programs,  "OPTEIGEN" , 
"DBEIGEN",  and  "PDPF” . Subroutines  in  the  software  package  "EISPACK" 
are  excluded  In  this  appendix  because  of  their  length.  Program  listings 
and  documentations  are  available  on  request. 


TrtlS  IS  PGM 
to  Hun  this 

Tu  HUN  This 
TJ  HUN  THIS 
This  IS  The 


' 06c  I-uEn. 

PG* : aRun  -lOao*na as:E  ISPaCh  un  mTs, 

huh;  bHUN  -L3AD+NAAS:tIS0AC^  UN  MfS 

PGM  In  uThER  uNI  V . SUoRuUTInE  PaChAgE  'EISPaCk'  IS  NEtDEO, 
NuMtHlCAL  ME  T rtOu  TU  GtT  EIGtN  VALuES  jF  MaTkIlES  03,Gu,ol 


PaThS, m, 13' ) 


"ns  is  i Ht  NuntHU.Au  Mt  i nuu  iu  ijt  i t u j t n vttuuto  jr  ■'lAiKiuto  Lio,bu»oi 
01  HENS  I JN  VARSClO),PHHlO),AllO),L(10J’,GOC10),Gl(10) 

REAL *6  ARClO,10),AlClu,10),wRllO),wiC10),ZRUU,lO),ZIClO,10), 

I SvlClO),Sv2t.lO),Sv3(iu) 

CUMPLEX  A,C 

STAHT  TO  GENERAT  MATRIX  Do  FOR  OOuduE  8EAMFuRMERS.»SSSSSbS.kiii.bS>biii3.S 
wRITECo, 104) 

l 0 4 FORMAT  C 3 X » • T T PE  IN  NO.  OF  ELEMENTS,  13') 

READC5, 105)  « 

105  FORMAUIS) 

hRITECo, 223) 

123  FORMATC3X, • TYPE  IN  NO.  OF  PaThS, m,  IS') 

HE AO ( 5 » 1 05 ) M 
RRITE(S*101) 

L 0 1 FORMA T( IX,  * TYPE  In  Ph I ( 1 ) / . . . PH  1 ( n ) , F l 4 . 0 ' ) 

REAO(b«10O)  CPHI  CKKM)  ,Ki\M=l,M) 

"HI  TE(8,  102) 

HEAD  15, 100)  CVASSCKKM) ,KKMsl ,M) , vAHN 
L 02  FORMAT! IX, • TYPE  In  V ARS C 1 ) , . . . v aRS ( M ) , V aHn , F 1 4 . 0 ' ) 

LOO  FORMATCF14.0) 

00  222  *m=1,m 

PHI CKM)sRhI (KM)*3.l4lb/l80. 

C(KM)sCMRLX(0.,3. 141b*C0S(PHlUM) ) ) 

A(KM)=CEXP!ClKM) ) 

222  CONTINUE 

m 0EN0TE  THE  NO.  UF  PAThS,  A 0ENOTE  ThE  nu.  UF  ELEMENTS 
00  5 11*1, K 
00  6 12=1 » K 
ARCI1, 12) *0. 

AICIt,I2)  =0. 
a CUNTlNUE 

3 COnTINuE 

DU  10  1*1, m 
OU  20  11*1, a 
00  21  Jl*l,K 

ARC  1 1 , J1  )sREAL(ACl  )<**(J1-I1))*VARSCI)+ARCI1,J1) 
aICI1,J1)*AIhaGCA(I)**CJ1-I1))*vAhS(I)>AICI1,J1) 
i 1 CONTINUE 

20  CUNTlNUE 

10  CONTINUE 

00  19  1=1, K 

19  "RITE (6, 117)  CARCI,J),Jst,K) 

00  29  1*1, k 

29  WRITECO, 1 17)  ( A I C I , J ) , J*  1 , * ) 

117  FORMAT(3x, /3020.S) 

Eno  JF  GENERATING  MATRIX  Oo  A************************************* 

RE aL ( 00 ) rs  AR;  AIMAGlDd)  IS  Al 

CALL  CGC10,K,AH,Al,"H,"I,0,ZR,ZI,Svl,Sv2,Sv3,lER) 

IF  (IER  .NE.  0)  hRITE  Co, 179) 

179  F0RMATC3X, 'ERROR  In  SU8  CG’) 

WRITE  Cts, Rod) 

*80  FORMATC3X, 'EIGENVALUES  OF  buUdLt  dE aMFURMER ' ) 
wRl  TE (b, 1 i 3 ) C*RCI),*ICI), 1*1, R) 

113  FuRMAT  CdX,2U20.8)  I^cXlCA^ 

OU  90  L*l»t\  _.,.e  wins  IS  V _ TMnfl 


iwb  — r > ' < * • » 

FuRMAT  CdX,2U20.8) 


OU  90  L*1,N 
uO(L)=vArn*hRCl) 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

Gl  (l)  = VARn*,»HiL)  *«RIL)**R(L)  iTJRttLSHED  TO  LDQ  

90  CUNT  IN  JE 

*RlTE(6,ll5)  (Gl  (Ll),Ll*l,K ) 

Us  FORmaRSa,  'EIGENVALUES  OF  Gl  ' * 3E2Q ) 

*RITE(6,U4)  (G0(Lt ),Ll*l,K) 

114  FURMARiX, 'EIGENVALUES  QF  Gu',3E2u.d) 

GO  TO  =»9 

999  KRITE(g#991) 

991  FURMA r (3X, 'METHOD  FAIL') 

99  STOP 

END 


C THIS  IS  PGM  'OPTEIGEN' 

C TO  RUN  THIS  PGM;.  *RuN  -lOAO*NaAS:E  ISPACa  Un  hTS. 

C TO  RUN  THIS  PGM  In  OThER  UNIV.  SUsRuUTInE  PACKAGE  'EISPaCK'  IS  NEEDED. 
C This  IS  The  NUMERICAL  WAT  TU  GET  THE  EIGENVALUES  JF  01  AND  G0,G1 
DIMENSION  VARSC10),PHI(10),AC1Q),C(10),GOCIu),G1(10) 

REALMS  ARUO,10),AI(10,10),*R(10),*I(10),Zh(10,10),ZI(10,1u) 

1 ,SVI(10),SV2(10),3V3C10) 

complex  A,C 

C START  TO  GENERATING  MATRIX  Q0-1Q1.  ********************************* 
*RITE(6, 104) 

104  FORMATC3A* ' TYPE  IN  NO.  OF  ELEMENTS  K,l3‘) 

RE AO  C 5* 105)  K 

105  FORMAT (13) 

*RIT£(6,223) 

223  FORMAT ( 3X , • TYPE  IN  NO,  OF  PATHS, 13') 

READ<5,105)  M 
WRITE16#  101  ) 

101  FORMA  R 1 X, ' TYPE  IN  PHI 1 1 ) , PHI  12 ) , . . . , PHI ( M ) , F 1 U . 0 ' ) 

READ  15  f 1 00 ) (PHI (KKM),KKM=1,M) 

*RI  TEC6,  102) 

READ ( 5 , 1 00  ) (VARS(KKM) ,KKMsl ,M) , VARN 

102  F0RMAT(2X# 'TYPE  IN  VARS ( 1 ),..., VaHS C M ), YARN, F 1 4 . 0 ' ) 

100  FORMAT ( F 1 4 , 0 ) 

DO  222  KM*l#M 

PHI  (KM)  sPhR KM)  *3,1416/180. 

C(KM)sCMPLX(0,#3. 1416*C0S(PHI (Km) ) ) 

A(KH)sCEXPCCCKM) ) 

222  CONTINUE 

C •*  DENOTE  TH£  NO,  OF  PATHS,  K DENOTE  THE  NO,  OF  ELEMENTS 
DO  5 I 1*1 ,K 
00  6 12*1, K 

IF  (II  .£0,  12)  ARCI1, I2)*VARN 

a I ( 1 1 , 1 2 ) *0. 

• CONTINUE 

CONTINUE 
Ow  10  1*1, n 
)u  29  Il«l,K 
. 21 

*-  ; . . i )**t*L(A(  I ) •*(  Jl-U  ) ) *VARS(  I ) tAR(  11 , Ji) 

. «*1-AG( A( I ) *•{ J 1-1 1 ) ) *VAKS( I ) +AI  ( II  , J 1 ) 

m * 1 •%»€ 

4 <•'  T iMt 
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11/  FQRMaT(/,3x,3020.8) 

C EnO  UK  GENERATING  QO“lJi  5 .&  .3  a 3 j>  i±  $ a 3 j>  3 a $ 4 i s S.i  a *3  a 3 a 5 a A & S a i a 5 3.5  a s a 
C REAL(Ql)  13  AR ; AIMaG(!)1)  IS  a! 

CALL  CG(l0,«,Art,Al,WH,«I,0,ZR,ZI,Svl,Sv2,Sv3,IER) 
lii  FORMAT (3x,2D£0.8) 

wRlTECo,  1 1 1 ) IEW 

111  FORM A I ( 2X » ' IE Rs  ' , I b ) 
wRITE(6*  1 1 2 ) 

112  fo«maT(2x, 'eigenvalues  of  u i ’ ) 

wRITECo*  113)  CwR(I),wI(I),I=l,iO 

C CALCULATE  EIGENVALUES  of  GO  and  Gl 
DO  90  L s 1 # * 

G0(L)*1 ,-VARn/wR(L) 

Gl  1L)s<nR(L)/VAhn-1  . 

90  CONTINUE 

wRITE(6,ll5)  (GlCLlJ.Llsl.KJ 

114  FORMAT  CSX, • EIGENVALUES  UF  GO  1 , 5E20 . o ) 

115  FORMaTCSx,  'EIGENVALUES  OF  G17',3E2o.o) 

wR I TE ( o * 1 14)  (G01L1),L1=1,K) 

STUR 

END 

C This  IS  PGM  * PDPF • . INPUT:  EIGENVALUES  of  G0.G1;  JUTPUT:PO.HF 
DIMENSION  ELTA(10),PQ(10},PF(lo) 

DO  99  IIIsl,20 
wRITECb.  103) 

READ  (5,100)  (ELTA(I),I*1#10) 

*RITE(6, 104) 

103  FORMAT (2X, ' TYPE  IN  10  THRESHOLD  VALUES') 

READCbilOo)  Gli,Gl2,G01»Go2 

100  FORMAT (F 1 4 . 0 ) 

104  FORMAT ( 2X , ' T YPE  IN  G 1 1 , G 1 2 , GO  1 » GO 2 , F14.u') 

CALL  PDPFCELTA/Gl l/G12/G0l,G02/PD/PF ) 

*HITE(6,200  ) Gil, G12,  GO  W GO  2 

WRITE (6, 30 l ) 

301  FURMATC2X, ' THIS  IS  PD  & PF ' ) 

WRITE (6,300)  lPD(J),PF(J),Jsi, 10) 

200  FURMAT(2X,2F18.6) 

300  FURmAT(2X,2F10.5) 

wRITE (6, 111) 

111  FQHMATC3A, 'STOP?  Y;  TYPE  NONZERO  INTEGER;  v:  T r PE  O') 

RE AO (5, 112)  ISTOP 

IF ( I S TOP  .NE,  0)  Gu  TO  199 

112  FURmA  T ( 1 3 ) 

49  CONTINUE 

1 99  STOP 


SUBROUTINE  PUPF(ELTA,Gll ,Gl2,G01,Go2,PD,PF  ) 

DIMENSION  ELTA(10),PO(lO),PF(10) 

DO  10  Ks  1 , 1 o 

P0U)*G1 1*EXP(-ELTA(K)/G1  l)/(Gll-G12)fGl2*ExP(-fcLTA(*)/ 
G12J/CG12-G1 1 ) 

PF (K )sG0 1 *E  X P ( “EL  T A ( K ) XG0l)/(G01-G02)+G02*ExP(-EL!A(K)/ 
G02)/(G02-G0l ) 

CONTINUE 
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Fig.  1.2  Ray  diagram  of  transmission  in  the  deep 
sound  channel  for  a source  on  the  axis 
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Figure  2.4  The  probability  density  function  of  phase  angle  as  a function  of  si 
nal-to-noise  ratio  in  dB 


Figure  2.5  90%  confidence  limits  on  measured  phase  as  a function  of  signal-to-noisa 

ratio  ir.  dB  • . 
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Figure  2.1  Computer-programmed  Phase-coherent  Demodulation  of 
CW  Transmission,  1550  h 3 Feb.  1965.  (Figure  16 
in  Reference  [4]) 
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Tj=  d^/c  : Time  delay  of  signal  from  first  path 
Ty*  d3/c  : Time  delay  of  signal  from  third  path 
where  dy=  D cos 

d3»  D cos  <J>3 
c:  sound  velocity 


Figure  3.1  Multipath  Signals  Geometry 
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(c)  Signal  estimate  of  1st  path  signal  (d)  Signal  estimate  of  2nd  path  signal 


Figure  3. 


3 (a)  Noise-free  Data  Received  at  1st  Element. 

(b)  Noise- free  Data  Received  at  2nd  Element. 

(c)  (d)  Phasor  Diagrams  of  Signal  Estimates  of  1st  and  2nd  path 
Signals. 
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Figure  5.3  Comparison  of  Optimal  performance  for  Detecting 
Two-path  Signals  with  Unequal  Strength. 

Signals  are  from  + 13 0 of  Horizontal  Axis. 


Figure  5*7  Comparison  of  the  Optimal  Performances  of  One,  Two,  and. 

Three-element  Arrays  for  Detecting  Two-path  Signals 
txon  * 15°  of  Horizontal  Axis.  Slgnal-to-nolse  Ratio 
Per  Element  Is  2 or  10. 


Figure  5.8  Performance  Comparison  Between  the  K-element  Array 

and  a Single  Phone  with  K Times  Signal- to-noise  Ratio 
Per  Element  of  Array  Case.  Signal- to-noise  Ratio  Per 
Element  (SNR)  Is  1. 


Figure  5*11  Performance  Comparison  Between  Optimal  Array  Detectors 

and  Single  Beamformer.  Two-path  Signals  Are  from  ± 15° 
of  Horizontal  Axis.  The  Signal-to-noise  Ratio  Per  Element 
(SNR)  Is  10. 


Figure  5*14  Performance  Comparison  Between  the  Optimal  Array  Detectors 
and  the  Double  Beamformer  Detectors*  Two -path  Signals 
Are  from  ± 15°  of  Horizontal  Axis.  The  Signal-to-noise 
Ratio  Per  Element  (SNR)  is  10. 
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